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Abstract 

By a generalized inverse of a given matrix, we mean a matrix that 
exists for a larger class of matrices than the nonsingular matrices, that 
has some of the properties of the usual inverse, and that agrees with 
inverse when given matrix happens to be nonsingular. In theory, there 
are many different generalized inverses that exist. We shall consider the 
Moore Penrose, weighted Moore-Penrose, Drazin and weighted Drazin 
inverses. 

New determinantal representations of these generalized inverse based 
on their limit representations are introduced in this paper. Applica¬ 
tion of this new method allows us to obtain analogues classical adjoint 
matrix. Determinantal representations for projection matrices A''’A, 
AA+, and A are presented as well. Using the obtained analogues of 
the adjoint matrix, we get Cramer’s rules for the least squares solution 
with the minimum norm and for the Drazin inverse solution of singular 
linear systems. The least squares solutions with the minimum norm 
and the Drazin inverse solutions of the matrix equations AX = D, 
XB = D and AXB = D are considered. Thus A B can be singular 
matrices of appropriate size. We get analogs of the Cramer rule for 
the generalized inverse solutions of these matrix equations and using 
their for determinantal representations of solutions of some differential 
matrix equations, X' -P AX = B and X' -P XA = B, where the matrix 
A is singular. 

Keywords: generalized inverse, Moore Penrose inverse, weighted 
Moore-Penrose inverse, Drazin inverse, weighted Drazin inverse, sys¬ 
tem of linear equation, matrix equation, differential matrix equation, 
Cramer’s rule. 
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1 Preface 


It’s well-known in linear algebra, an n-by-n square matrix A is called in¬ 
vertible (also nonsingular or nondegenerate) if there exists an n-by-n square 
matrix X such that 

AX = XA = In. 

If this is the case, then the matrix X is uniquely determined by A and is 
called the inverse of A, denoted by A“^. 

By a generalized inverse of a given matrix, we mean a matrix that exists 
for a larger class of matrices than the nonsingular matrices, that has some of 
the properties of the usual inverse, and that agrees with inverse when given 
matrix happens to be nonsingular. 

For any matrix A € consider the following equations in X: 


AXA = A; 

XAX = X; 
(AX)* = AX; 

(XA)* = XA. 

and if m = n, also 

AX = XA; 

A^+^X = A^ 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


For a sequence Q of {I, 2,3,4,5} the set of matrices obeying the equations 
represented in Q is denoted by A{Q}. A matrix from A{0} is called an 
^-inverse of A and denoted by A^®^ 

Consider some principal cases. 

If X satisfies the all equations ([I])-(j3]) is said to be the Moore-Penrose 
inverse of A and denote A+ = MoorePenrose inverse was 

independently described by E. H. Moore [I] in 1920, Arne Bjerhammar [2] in 
1951 and Roger Penrose [3] in 1955. R. Penrose introduced the characteristic 
equations ([I])-(j4]). If det A ^ 0, then A+ = A“^. 

The group inverse A^ is the unique Ad>2>5) inverse of A, and exists 
if and only if IndA = min{A: : rankA^"*"^ = rank A^} = 1. 

A matrix X = A^ is said to be the Drazin inverse of A if Q (for 
some positive integer k), ([2]) and ([5]) are satisfied, where k = IndA. It is 
named after Michael P. Drazin [3]. In particular, when IndA = 1, then the 
matrix X is the group inverse, X = A®. If IndA = 0, then A is nonsingular, 
and A^ = A~^. 
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Let Hermitian positive definite matrices M and N of order m and n, 
respectively, be given. For any matrix A € the weighted Moore- 

Penrose inverse of A is the unique solution X = A^ ^ of the matrix 
equations ([I|) and ([2]) and the following equations in X [5]: 

(3M) (MAX)* = MAX; (4A) (NXA)* = NXA. 


In particular, when M = 1^ and N = I„, the matrix X satisfying the 
equations ([I|), (l2|), (3M), (4N) is the Moore-Penrose inverse A+. 

The weighted Drazin inverse is being considered as well. 

To determine the inverse and to give its analytic solution, we calculate 
a matrix of cofactors, known as an adjugate matrix or a classical adjoint 
matrix. The classical adjoint of A, denote Adj[A], is the transpose of the 
cofactor matrix, then A~^ = —Representation an inverse matrix by 
its classical adjoint matrix also plays a key role for Cramer’s rule of systems 
of linear equations or matrices equations. 

Obviously, the important question is the following: what are the ana¬ 
logues for the adjoint matrix of generalized inverses and, consequently, for 
Cramer’s rule of generalized inverse solutions of matrix equations? 

This is the main goal of the paper. In the paper we shall adopt the 
following notation. Let be the set of m by n matrices with complex 

entries, be a subset of in which any matrix has rank r, 1^ be 

the identity matrix of order m, and ||.|| be the Frobenius norm of a matrix. 

Denote by a j and aj. the jth column and the zth row of A G 
respectively. Then a*- and a* denote the jth column and the zth row of a 
conjugate and transpose matrix A* as well. Let (b) denote the matrix 
obtained from A by replacing its jih. column with the vector b, and by 
Aj. (b) denote the matrix obtained from A by replacing its ith row with b. 

Let a := {ai,...,ak} C {!,..., m} and /3 := {/3i,...,/3fc} C n} 

be subsets of the order 1 < A: < min {m,n}. Then A^ denotes the minor 
of A determined by the rows indexed by a and the columns indexed by (3. 
Clearly, |A“| denotes a principal minor determined by the rows and columns 
indexed by a. The cofactor of aij in A € is denoted by g|-|A|. 

For 1 < A < n, := { a : a = (ai,..., a^), 1 < ai < ... < < n} 

denotes the collection of strictly increasing sequences of k integers chosen 


from {1,..., n}. 
1 < p < A, let 


Let A/j :— Lk,m ^ ^k,n- For fixed o. G Lp^rn-, /3 G Lp,m 

Ik,m {cn) ■— {L . I G Lk, mil 2 Q;}, 

Jk,n (/3) := {J : J G Lfc^n, J D /?}, 

-Afc (q:,/3) .— Ik,m {oi) ^ Jk,n (/3) 
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For case i G a and j & (3, we denote 

Ik,m{^} ■— {ck . Ct G G Cv}, .— {/3 . /3 G £ /3}) 

The paper is organized as follows. In Section 2 determinantal representations 
by analogues of the classical adjoint matrix for the Moore Penrose, weighted 
Moore-Penrose, Drazin and weighted Drazin inverses are obtained. 

In Section 3 we show that the obtained analogues of the adjoint matrix 
for the generalized inverse matrices enable us to obtain natural analogues of 
Cramer’s rule for generalized inverse solutions of systems of linear equations 
and demonstrate it in two examples. 

In Section 4, we obtain analogs of the Cramer rule for generalized inverse 
solutions of the matrix equations, AX = B, XA = B and AXB = D, 
namely for the minimum norm least squares solutions and Drazin inverse 
solutions. We show numerical examples to illustrate the main results as well. 

In Section 5, we use the determinantal representations of Drazin inverse 
solution to solutions of the following differential matrix equations, X' + 
AX = B and X' + XA = B, where A is singular. It is demonstrated in the 
example. 

Facts set forth in Sections 2 and 3 were partly published in [6] , in Section 
4 were published in [7l[8] and in Sections 5 were published in [8]. 

Note that we obtained some of the submitted results for matrices over 
the quaternion skew field within the framework of the theory of the column 
and row determinants ( IM]). 

2 Analogues of the classical adjoint matrix for the 
generalized inverse matrices 

For determinantal representations of the generalized inverse matrices as ana¬ 
logues of the classical adjoint matrix, we apply the method, which consists on 
the limit representation of the generalized inverse matrices, lemmas on rank 
of some matrices and on characteristic polynomial. We used this method 

at first in [6] and then in [8]. Liu et al. in m deduce the new determi- 

( 2 ) 

nantal representations of the outer inverse A^^^ based on these principles 
as well. In this chapter we obtain detailed determinantal representations by 
analogues of the classical adjoint matrix for the Moore Penrose, weighted 
Moore-Penrose, Drazin and weighted Drazin inverses. 
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2.1 Analogues of the classical adjoint matrix for the Moore 
- Penrose inverse 


Determinantal representation of the Moore - Penrose inverse was studied 
in m, ffMn. The main result consists in the following theorem. 

Theorem 2.1 The Moore - Penrose inverse = (a^) € of A £ 

(j^mxn following determinantal representation 



E 

{a,/3)eNr{j,i} 




d 

Odj i 



E 

( 7 , s)eNr 




^ ^ i, j < n. 


This determinantal representation of the Moore - Penrose inverse is based 
on corresponding full-rank representation m- ifA = PQ, where P € 
and Q G then 

A+ = Q*(P*AQ*)“^PT 

For a better understanding of the structure of the Moore - Penrose inverse 
we consider it by singular value decomposition of A. Let 

AA*Uj = afui, i = l,m 
A*Avi = afwi, i = l,n, 

and the singular value decomposition (SVD) of A is A = USIV*, where 


U= [uiU 2 ...u^] €C™x-, U*U = I™, 
V= [V1V2...V„] GC-X-, V*V = I„, 


5] = diag(cTi,U2,...,u^) 

Then [3], A+ = VSl+U*, where = diag((Tj“^, 

We need the following limit representation of the Moore-Penrose inverse. 


Lemma 2.1 /gg|/ If A £ then 

A+ = lim A* (AA* + AI)"^ = lim (A*A + AI)”^ A*, 

A->0 A-^0 

where A G M+, and M+ is the set of positive real numbers. 

Corollary 2.1 J2^ If A £ then the following statements are true. 
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i) //rank A = n, then A+ = (A*A)~^ A* . 
n) //rank A = m, then A+ = A* (AA*)“^. 
in) /frankA = n = m, then A"*" = A~^ . 

We need the following well-known theorem about the characteristic poly¬ 
nomial and lemmas on rank of some matrices. 

Theorem 2.2 124^ Let dr he the sum of principal minors of order r of 

A € Then its characteristic polynomial pA (t) can be expressed as 

Pa {t) = det {tl — A) = t^ — dit'^~^ + d 2 t^~^ — ... (—1)"' dn- 

Lemma 2.2 //As then rank (A* A) ■ < r. 

Proof. Let Pj^ {—Ojk) S (k i), be the matrix with —Ojk in the 

(f, k) entry, 1 in all diagonal entries, and 0 in others. It is the matrix of an 


elementary transformation. It follows that 



(A*A).,(a^)-nP*'«(-«i'«) = 

( E ■ 

k¥=j 

■ "t ■ 

El 

k¥=j 




k^i 

E “nfcafcl ■ 
V k^j 

■ <, ■ 

i—th 

E ankakn 

k^j / 


The obtained above matrix has the following factorization. 

/ E 


V 


an ai2 • • • «r 

0-21 ®22 • • • 


^nl nn2 • • • ®n 




aij ■ 

* 

a-nj ■ 

i—th 
( On 
0 

\ Q-ml 


El 

k¥=j 

E a^j^akn 

k¥=j 


0 

1 

0 

i—th 


0 

^mn / 


j - th. 


6 
















Denote by A := 


/ On 


V Oml 


0 

1 

0 

i—th 


^In ^ 




j — th. The matrix A is ob¬ 


tained from A by replacing all entries of the jth row and of the ith column 
with zeroes except that the (j, i) entry equals 1. Elementary transforma¬ 
tions of a matrix do not change its rank. It follows that rank (A* A) 




min I rank A*, rank A |. Since rank A > rank A = rank A* and rank A* A = 
rank A the proof is completed.■ 

The following lemma can be proved in the same way. 


Lemma 2.3 If A £ then rank (AA* 


:>■/ — 


< r. 


Analogues of the characteristic polynomial are considered in the following 
two lemmas. 


Lemma 2.4 If A £17' 


and A G M, then 


det((AI„ + A*A) Ja*.)) 


= + ... + 


.(b') 


(7) 


where Cn'^^ = 


(A*A) - (a* ) and = Y. 


/3eJs.n{i} 


A*A),, a^ 


for 


all s = l,n — 1, i = 1, n, and j = 1, m. 


Proof. Denote A*A = V = {vij) £ Consider (AI^ -|- V) ■ (v,j) G 

Cnx”’. Taking into account Theorem 12.21 we obtain 


|(AI„ + V) . (v,,)| = diA”-i + d 2 \^-^ + ... + dn, (8) 

where dg = Y l(A)^| is the sum of all principal minors of order s 
0eJs,n{i} 

that contain the i-th column for all s = l,re — 1 and d^ = detV. Since 

V. j = Y where a* is the Th column-vector of A* for all I = l,n, then 

i 

we have on the one hand 


|(AI + V) , (v.,)| = E l(AI + V)., (a>,,)| = 
EI(AI + V),aa*,)|-a,, 
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Having changed the order of summation, we obtain on the other hand for 
all s = 1, n — 1 


4= E (V)^ = E E(V.(a>H))'’ 

/3GJs.n{i} 0&Js,n{i} I 




E E (V..(a:,))^ 

I /SGJs, 7l{i} 


( 10 ) 




By substituting ([9]) and (fTOj) in ([8]), and equating factors at an when I = j, 
we obtain the equality ([7t). ■ 

By analogy can be proved the following lemma. 


Lemma 2.5 // A G 


and A G M, then 


det ((AI^ + AA*),-, (a*)) = + ... + 


where rm^ = |(AA*)j, (a* )| and = Y. i((AA*)j, (a* ))"| for all 

flG/s,7n{H 


Ab') 


s = 1, n — 1, i = l,m, and j = 1, n. 

The following theorem and remarks introduce the determinantal represen¬ 
tations of the Moore-Penrose by analogs of the classical adjoint matrix. 


Theorem 2.3 // A G and r < min{m, n}, then the Moore-Penrose 

inverse A+ = G possess the following determinantal represen¬ 


tations: 


or 


I ^^Jr, n'{^} 

% = - 


E (A* A).. 




| Q:G/r,r7i{j} 


E (A*A)^ 

n 


E |((AA*),-,(a*)) 


( 11 ) 


% = 


E l(AA* 

Ot^Ir, m 


( 12 ) 


for all i = l,n, j = 1, m. 


Proof . At first we shall obtain the representation (jllh . If A G K_|_, then 
the matrix (AI-|-A*A) G is Hermitian and rank (AI-1-A* A) = n. 

Hence, there exists its inverse 


(AI +A*A)~^ 



/ Lii 

L 2 I 

Ln 1 

\ 

1 

Li2 

L 22 

■ ■ Ln2 


det (AI -1- A*A) 






\ Ll n 

L2n ■ 

• • Lfifi 

/ 


8 





















where Lij {yi,j = l,ra) is a cofactor in AI + A*A. By Lemma l2.ll A"^ = 
lim (AI + A*A)~^ A*, so that 

A->-0 


A"*" = lim 
A^o 


/ det(AI+A*A) 

det(AI+A*A) 

det(AI+A*A) „(a*J 
\ det(AI+A*A) 


From Theorem 12.21 we get 


det(AI+A*A) 

det(AI+A*A) 


det(AI+A*A) „(a*„) 
det(AI+A*A) / 


det (AI + A*A) = A” + + ... + dn, 


(13) 


where (Vr = l,n — l)isa sum of principal minors of A* A of order r and 
dn = det A* A. Since rank A* A = rank A = r, then dn = dn-i = ... = 
dr+i = 0 and 

det (AI + A*A) = A^ + diA^-^ + + ... + (14) 


In the same way, we have for arbitrary 1 < i < n and 1 < j < m from 
Lemma 12.41 


det (AI + A*A) . (a* ) = + ... + 


where for an arbitrary 1 < k < n — 1, ^ 

^^Jk, n{i} 

In^'^ = det (A*A) ■ (a* y By Lemma [221 rank (A*A) * (a* ) < r so that 


A*A).i(aL 




and 


if /c > r, then 


A*A),i(a* 


= 0, (V/3 G Jk,n{i}yi = l,ra,Vj = l,m). 


Therefore if r + 1 < fc < n, then ^ 

I^^Jk, ri{^} 

liij) 


A*A),iia*. 




= 0 and 


In^ = det (A*A) ^ (^a*y = 0, (Vi = l,n, Vj = l,m). Finally we obtain 

det (AI + A*A) . (aT) = + ... + (15) 

By replacing the denominators and the numerators of the fractions in 
entries of matrix (1131) with the expressions (fT4)) and (fT5]l respectively, we get 


A”*^ = lim 
A^o 


A"+diA"-i+...+dr.A’*-^ 


4"i)A"-i+...+i("i)A"-’' 

A"+ciiA"-i+...+dr.A’»-'' 


4i"*)A"-i+...+i(.i"*)A"-’' \ 

A'*+(iiA"-i+...+drA"-'' 
A'*+(iiA"-i+...+dr.A"-'’ / 
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,( 11 ) 


,(nl) 


(Im) \ 

/ 


di'p 


From here it follows m- 

We can prove m in the same way. ■ 

Corollary 2.2 If A ^ and r < min {m,n} or r = m < n, then the 

projection matrix P = A^A can be represented as 


P = 


Pij 


dr (A*A) 


where d j denotes the jth column of (A*A) and, for arbitrary 1 < i,j < n, 
Pij = E ((A*A).i(d.j))^ ■ 

/3G i/r,n 

Proof. Representing the Moore - Penrose inverse A"*" by (|lip . we obtain 


1 


dr {A*A) 


/hi 

h2 ■ 

llm^ 


/ an 

012 

ain '' 

hi 

h2 ■ 

^2m 


021 

022 

02n 

\lnl 

ln2 

Inm / 


\Om 1 

am 2 

am nj 


Therefore, for arbitrary 1 < < n we get 

Pij = E E ((A*A).i(a*J)^ 

k /3G</r, n 

= E E = E 


/3€Jr, niO k 




^kj 


A*A).i(d* 


Using the representation (|12l) of the Moore - Penrose inverse the following 
corollary can be proved in the same way. 

Corollary 2.3 If A G where r < min {m, n} or r = n < m, then a 

projection matrix Q = AA+ can be represented as 


Q = 


Qij 


dr (AA* 


where g, denotes the ith row of (AA*) and, for arbitrary I < i,j < m, 

9i,= E |((AA-),.(g,));|. 


aG/r,m{i} 
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Remark 2.1 //rank A = n, then from Corollarv \2. l\ we aet A.^ = (A* A) ^ A*. 
Representing (A*A)“^ by the classical adjoint matrix, we have 

/det(A*A).i(a*i) ... det(A*A),i (a*^) 

A+ = 

det(A*A) \^det(A*A),„(a*^) ... det(A*A),„ (a*^) 

If n < m, then m is valid. 

Remark 2.2 As above, i/rank A = m, then 



A+ = 


'det(AA*)i.(at ) 


det(AA-) 


If n> m, then r/gj) is valid as well. 


det(AA*)^,(a^ )\ 

• ( 17 ) 

det(AA*)^, (a; )/ 


Remark 2.3 By definition of the classical adjoint Adj{A) for an arbitrary 
invertible matrix A G Qng f^ay put, Adj{A) ■ A = det A • R. 

// A G and rank A = n, then by Corollary \2.1[ A+A = In- 

Representing the matrix A"*" by as A+ = det(A*A) ^ obtain LA = 

det (A*A) • In- This means that the matrix L = {hj) G is a left 

analogue of Adj{A) , where A G and kj = det(A*A).i for all 


i = l,n, j = 1, m. 

//rank A = m, then by Corollary \2.1i AA+ = 1^. Representing the 
matrix A~^ by JT^ as A"*" = det(AA*) ^ obtain AR = Im ■ det(AA*). 
This means that the matrix H = (rij) G is a right analogue ofAdj{A), 

where A G and rij = det(AA*)j, (a/) for all i = l,n, j = l,m. 

If A £ and r < m.m{m,n}, then by 1771) we have A^ = ^ (a*a) ^ 


where L = (/jj) G C' 


and lij 


((A*A) ,(a*.))^| for all 


= E 

13^ Jr, n{i} 

i = l,n, J = l,m. From Corollary Ig.gl we get LA = dr (A*A) ■ P. The 
matrix P is idempotent. All eigenvalues of an idempotent matrix chose from 
1 or 0 only. Thus, there exists an unitary matrix U such that 


LA = dr {A*A) Udiag (1,..., 1,0,..., 0)U*, 

where diag (1,..., 1,0,..., 0) G is a diagonal matrix. Therefore, the 

matrix L can be considered as a left analogue of Adj{A), where A G 

In the same way, if A £ and r < mm{m,n}, then by m we have 

7^^ = dr {A A*) ’ ^ ^ Vij = Y. l((AA*)j. (a* )) “I/or 

a&Ir,m{j} 
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all i = l,n, j = l,m. From Corollary \ 2. 3\ we get AR = dr (AA*) • Q. The 
matrix Q is idempotent. There exists an unitary matrix V such that 

AR = dr (AA*) Vdiag (1,..., 1, 0,... , 0)V*, 

where diag (1,..., 1,0,... , 0) € Therefore, the matrix Yl can be 

considered as a right analogue of Adj{ A) in this case. 

Remark 2.4 To obtain an entry of A^ by Theorem \2.1\ one calculates 
{CJ^Cfn + determinants of order r. Whereas by the equation 

m we calculate as much as (C^ + determinants of order r or we 

calculate the total of (C^ + determinants by /fil)) . Therefore the 

calculation of entries of A"*" by Theorem \2.3\ is easier than by Theorem \2.1\ . 


2.2 Analogues of the classical adjoint matrix for the weighted 
Moore-Penrose inverse 


Let Hermitian positive definite matrices M and N of order m and n, re¬ 
spectively, be given. The weighted Moore-Penrose inverse X = A)[^ ^ can 
be explicitly expressed from the weighted singular value decomposition due 
to Van Loan [25] 

Lemma 2.6 Let A € There exists U G y g £^nxn gQfj^gfying 

U*MU = and V*N-iV = R such that 


A = U 


D 0 
0 0 


V* 


Then the weighted Moore-Penrose inverse AJ^^ can be represented 

At,A. = N-‘V ( “■' o)u*M, 

where D = diag{ai,a 2 , ...,ar), cri > a 2 > ... > ar > 0 and is the nonzero 
eigenvalues o/N“^A*MA. 

For the weighted Moore-Penrose inverse X = A)[^^ we have the following 
limit representation 

Lemma 2.7 ( J2S\/ . Corollary S.f.) Let A G A^* = N~^A*M. Then 

A+,^ = lim(AI + A»A)-iA# 
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Denote by a ^ and aj the jth column and the ith row of respectively. 
By putting instead A*, we obtain the proofs of the following two lemmas 
and theorem similar to the proofs of Lemmas 12.21 12.41 and Theorem 12.31 
respectively. 

Lemma 2.8 If A'^A G then 


rank 




< rank 



for all i = l,n and j = I, m 


Analogues of the characteristic polynomial are considered in the follow¬ 
ing lemma. 


Lemma 2.9 // A G (^"ix^ and A G K, then 
det ^ + A**A^ 




where Cn^'^ = 


[A#A)^ 


1 a, 

I \ -j 


and 

pG.Js, 


({A«A),.(a«,)) 


for 


all s = l,n — 1, i = 1, n, and j = 1, m. 


The following theorem introduce the determinantal representations of 
the weighted Moore-Penrose by analogs of the classical adjoint matrix. 


Theorem 2.4 //A G and r < min{m,n}, then the weighted Moore- 

Penrose inverse A^^ = possess the following determinantal 

representation: 


% = 


E , 


A#A) . 
/ . 1 




E (A«A) 


(18) 


^^Jr, n 


for all i = l,n, j = 1, m. 


2.3 Analogues of the classical adjoint matrix for the Drazin 
inverse 

The Drazin inverse can be represented explicitly by the Jordan canonical 
form as follows. 
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Theorem 2.5 


|^//A G C"""” 
A = P 


with IndA = k and 




p-i 


where C is nonsingular and rank C 
k, then 



= rank A^, and N is nilpotent of order 



(19) 


Stanimirovic’ |28] introduced a determinantal representation of the Drazin 
inverse by the following theorem. 

Theorem 2.6 The Drazin inverse A^ = of an arbitrary matrix A G 

(j^nxn possesses the following determinantal representation 


E 

{a,g)eNr^.{j,i} 


(A" 


(tdj i 


Oij 


E (A^ 


(7, <5) 67V. 


I 


1 < hj < n\ 


( 20 ) 


where s > k and r^. = rankA^. 


This determinantal representations of the Drazin inverse is based on a full- 
rank representation. 

We use the following limit representation of the Drazin inverse. 
Lemma 2.10 If A £ then 

A^ = lim (Xln + A^+^'l A^, 

A^O V / 

where k = IndA, A G M+, and M+ is a set of the real positive numbers. 
Since the equation Q can be replaced by follows 

(6a) XA^+i = A^ 


the following lemma can be obtained by analogy to Lemma 12.101 
Lemma 2.11 If A G then 

A^ = lim A^ (Xln + A^+^'l , 

A^-O V / 

where k = IndA, X G M+, and M_|_ is a set of the real positive numbers. 
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(k) (k) h 

Denote by ' and ' the jth column and the ith row of A respectively. 
We consider the following auxiliary lemma. 

Lemma 2.12 If A € with IndA = k, then for all i,j = l,n 

rankAf+i < rankA^+^ 

Proof. The matrix may by represent as follows 

/A (fc) A (fc) \ 

/ A «is«k • • • A aisalri \ 

S=1 S=1 


A (fc) A (k) 

\ CLnsO^sl • • • A (^ns^sn j 

\ S=1 S=1 / 

Let Pii{—aij) G (Z 7^ i), be a matrix with —aij in the (Z,i) entry, 

1 in all diagonal entries, and 0 in others. It is a matrix of an elementary 
transformation. It follows that 

/ A (fc) A (fc) \ 

/ 2 ^ aisalf ... 2 ^ aisok \ 




Y1 ^nsCl^sl 




The obtained above matrix has the following factorization. 


A 

s¥=j 


(Jiscis 


A (fc) 


A (fc) 

/ . OinsO^sn 


fail ... 

0 


0 ... 

1 

... 0 

\aril 

0 

• • • Ojyin/ 



/ 



Jk) 


Oil 

®12 

• ®ln 

Jk) 

Ak) 


®21 

®22 

• ®2n 



„(fc) 

\®nl 

®n2 
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Denote the first matrix by 


A := 


/ On 


\ Ofil 


0 

jth 


Oln 


Onn / 


ith. 


The matrix A is obtained from A by replacing all entries of the ith row 
and the jth column with zeroes except for 1 in the (i, j) entry. Elemen¬ 
tary transformations of a matrix do not change its rank. It follows that 
rankA^"*"^ ~ min |rankA^, rank a|. Since rank A > rankA^ the 

proof is completed. ■ 

The following lemma is proved similarly. 

Lemma 2.13 // A G yjith JndA = k, then for all i,j = l,n 

rank A^+i < rankA^+^ 


Lemma 2.14 If A £ C"'Xn A G M, then 

det ((AI„ + A^+i)j, (af ^)) = + ... + (21) 


where rn^'^ = 


and = Y 

(A‘+-(a!‘>))" 



aGL.nli} 



1 , n — 1 and i,j = 1, n. 


Proof. Consider the matrix ^(AI^ + 
account Theorem 12.21 we obtain 


Ain + A 




= diA”-^ + d2X^‘ 


G Taking into 

^ + dn, ( 22 ) 


where dg = is the sum of all principal minors of order 

“eT.nlil 

s that contain the j-th row for all s = l,n — 1 and dn = det A^"*"^. Since 
^(fc-i-i) _ Yfajia.^^\ where a^^ is the Ith row-vector of A^ for all I = l,n, 

then we have on the one hand 


(^(Aln + A*^+^)j.(aj ^)^ (AH- A^+^)^ 


Ea,r (AI + A^+i), a; 




(23) 
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Having changed the order of summation, we obtain on the other hand for 
all s = 1, n — 1 


4 = E E E 

a&Is,n{j} Ct&h,n{j} I 

Za.r E 

I Ol^Is,n{j} 






(24) 


By substituting (1251) and (1241) in (1221) . and equating factors at aji when I = i, 
we obtain the equality (1211) . ■ 


Theorem 2.7 If InclA = k and rank = rank = r < n for A G 
then the Drazin inverse A^ = possess the following 

determinantal representations: 




E 

aG/r,n{i} 



E KA-'+'e 

(AGlr,n 


a 

a 


and 


for all i,j = 1, n. 


hj 


E , 

/3 G t/r, Ti { ^ } 



E 

PeJr ,n 


(A‘+>)? 


d 

d 


(25) 


(26) 


Proof. At first we shall prove the equation (j25p . 

If A € M+, then rank (AI + A^+^) = n. Hence, there exists the inverse 
matrix 



/ I?ii 

R 2 I 

• • Rn 1 

\ 


Ri2 

R 22 

• • Rn 2 


V J det (AI + A^+i) 






\ Rln 

R 2 n 

■ ■ Rn n 

/ 


where Rij is a cofactor in AI + A*^+^ for all i,j = l,n. By Theorem 12.111 

A^ = lim A^ (AI„ + A^"''^) so that 
A^-o ^ ' 


A^ 


A™odet (AI + A^+i) 




(k) p \ 
Z^s=l “is ^ns 

\^n (k) p 

/_^s=l rins I 
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det(AI+A'=+l)^ det(AI+A''+l)^ \ 

det(AI+A'=+i) ■ ■ ■ det(AI+Afc+i) 

det(AI+A'=+l)^ det(AI+A''+l)^ 

det(AI+A''+i) ■ ■ ■ det(AI+A''+i) / 

Taking into account Theorem 12.21 , we have 

det (^AI + = A” + diA^-i + daA'^"^ + ... + dn, 


lim 

A^o 


(27) 


where dg = Y1 I ^ principal minors of of 

OL^Is,71 

order s, for all s = l,n — 1, and dn = det A^"*"^. Since rankA^"*“^ = r, then 
dn = dn-i = ... = dr+i = 0 and 

det (aI + A^+^) = A" + diA”-^ + d 2 X^-^ + ... + drA”“^ (28) 


By Lemma 12.141 for all i,j = l,n, 

det (aI + A^+i) + ... + 


where for all s = 1, n — 1, 



0 'Gls,n{j} 


and = detAj+i 


By Lemma [2.121 rankA^"*"^ ) < r, so that if s > r, then for all 
« G Is,n{i} and for all i,j = l,n, 


k+l f(k) 


3 ■ 


= 0 . 


Therefore if r + 1 < s < n, then for all i^j = 1, n, 


i?''= E 

“G/s.n{i} 




= 0 , 


and In^'^ = det A^^^ obtain 

det (^AI + A*^+^) (^af + ... + (29) 
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By replacing the denominators and the nominators of the fractions in the 
entries of the matrix (j27p with the expressions ()28p and (|29l) respectively, 
finally we obtain 




lim 

A->-0 


A'*+diA"-i+...+drA"-’' 

\ A’*+diA"-i+...+(irA"-^ 


;,(!")A"-i+...+«^i")A'*-’- \ 

A"+diA'*-i+...+(irA"-’' 
A"+(iiA'*-i + ...+(irA"-’' / 



where for all = l,n, 



E 

«e^r,n{i} 


A fc +1 
^3- 


a; 


(fc) 


, dr 


E !(a-): 


The equation (1261) can be proved similarly. 

This completes the proof. ■ 

Using Theorem 12.71 we evidently can obtain determinantal representa¬ 
tions of the group inverse and the following determinantal representation of 
the identities A^A and AA^ on R{A^) 

Corollary 2.4 If IndA = 1 and rankA^ = rank A = r < n for A G 
then the group inverse A^ = £ £^nxn pggg^gg following 

determinantal representations: 


for all i,j = 1, n. 


E 


n ^£U,Ti{t} 

<3 = - 


Af. (a. 


E l(A2)^ 

oe/r.n 


= 

% 


E 

/3GJr,n{j} 


(At (aj));. 


E (At? 

/3GJr ,n 


(30) 
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Corollary 2.5 If IndA = k and rankA^^^ = rankA^ = r < n for 
A G then the matrix AA^ = {qij) € possess the following 

determinantal representation 


OieIr,n{j} 

qij = 

<AGlr,n 

for all i,j = 1, n. 

Corollary 2.6 If IndA = k and rankA^+^ = rankA^ = r < n for 
A G then the matrix A^A = (pij) G possess the following 

determinantal representation 






for all i,j = 1 , n. 

2.4 Analogues of the classical adjoint matrix for the W- 
weighted Drazin inverse 

Cline and Greville [30] extended the Drazin inverse of square matrix to 
rectangular matrix and called it as the weighted Drazin inverse (WDI). 
The W-weighted Drazin inverse of A G with respect to the matrix 

W G is defined to be the unique solution X G of the following 

three matrix equations: 

1) (AW)*^+iXW = (AW)^ 

2) XWAWX = X; (33) 

3) AWX = XWA, 

where k = max{/nd(AW),/n(i(WA)}. It is denoted by X = A^^w- In 
particular, when A G and W = 1^ , then A^ ^r reduce to A^ . If 

A G is non-singular square matrix and W = 1^, then Ind{A) = 0 

and Ad^w = = A“^. 

The properties of WDI can be found in (e.g., [3IHM1). We note the 
general algebraic structures of the W-weighted Drazin inverse [M]. Let for 
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A e and W € exist L € and Q € C"""" such that 


A = L 


An 0 
0 A 22 


Q ^ W = Q 


Wii 0 
0 W22 


L -1 


Then 


^d,W = L 


(WnAnWn)-i 0 \ 1 

0 0 y ^ ’ 


where L, L, An, Wn are non-singular matrices, and A 22 , W 22 are nilpotent 
matrices. By [29] we have the following limit representations of the W- 
weighted Drazin inverse. 


Ad,w = lim (aI™ + (AW)'=+2) ' (AW)'=A (34) 

and 

Ad,w = limA(WA)'^ (aI„ + (WA)'=+2) (35) 

where A G M+, and M+ is a set of the real positive numbers. 

Denote WA =: U and AW =: V. Denote by and the jth 

column and the ith. row of respectively. Denote by := (AW)*^A G 
^mxn ^ ^ WAW G 


Lemma 2.15 If AW = V = (vij) G with Ind'V = k, then 

rank . (36) 


Proof We have V'=+2 = V'^W. Let Pis{-Wjs) G (s / i), 

be a matrix with —Wjs m the {i,s) entry, 1 in all diagonal entries, and 0 
in others. The matrix Pj^ (—Wjs), {s 7 ^ i), is a matrix of an elementary 
transformation. It follows that 




_(fc) 

• ■ 

—{k) — 

■ L v\fWsm 



-{k) - 
VmiWsl . 

V 

'^raj 

-{k} - 

■ VinsWsm 



i—th 



We have the next factorization of the obtained matrix. 


( E ■ 

-ik) 

■ •• 

-{k) - 

. L ns '^sm 



-{k) - 

2_j VmsWsl . 

V 

rnj 

-{k) - 

^+3 


i—th 
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"11 


12 




"21 


"22 


v 




ml 


Denote W := 


m2 


/ 1(^11 


\ Wnl 


yik) \ 
'^In 
_(fc) 


_(fc) 

Vmn / 


0 

1 

0 

i—th 


( Wll 
0 

V Wnl 


0 

1 

0 

i—th 


'^Im 

0 

'^nm. 


j - th. 


Wlm 

0 

Wri 


j — th. The matrix W is ob- 


/ 


tained from W = WAW by replacing all entries of the jth row and the 
ith column with zeroes except for 1 in the {i,j) entry. Since elementary 
transformations of a matrix do not change a rank, then rank — 

min < rank V^, rank W >. It is obvious that 


rankV^ = rank(AW)^A > rank (AW)*^+^, 
rankW > rank WAW > rank (AW) ^+2. 


From this the inequality (l.ldjl follows immediately. ■ 

The next lemma is proved similarly. 

Lemma 2.16 //WA = U = (uij) G with IndJJ = k, then 

rank < rank (^U^+ 2 ) , 

where := A{WA)^ € C™""" 

Analogues of the characteristic polynomial are considered in the follow¬ 
ing two lemmas. 

Lemma 2.17 If AW = V = (vij) G u^ith InclY = k and A G M, 

then 


Xlm + V 


k+2 


_(fc) 

j 


= + ... + 


( 37 ) 


where Cm^ = det (V^+2) ^ andci^^^ = ^ det f(V^+2) ^ 

m{T 

for all s = l,m — 1, i = 1 , m, and j = 1 , n. 
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Proof. Consider the matrix (AI + S Taking into 

account Theorem E 2 ] we obtain 

(aI + V^+ 2 ) ^ I = diA™-! + (isA™-^ + ... + dm, (38) 

where dg = I is the sum of all principal minors of order 

d^Js. m{d 

s that contain the z-th column for all s = l,m — 1 and dm = det 

/ \ 


Since = 


V 


—{^) — 
llvliwii 
i 


-{k) - 


= where is the Zth column- 

l 


-Xk) 


vector of = (AW)^A and WAW = W = (wu) for all I = l,n, then we 
have on the one hand 


(AI + V^+2) ^ I = E I (AI + V^+2) ^ 


E (AI + V"+2) 


-fk) 


(39) 


Wli 


Having changed the order of summation, we obtain on the other hand for 
all s = 1 , m — 1 


4 = E (v‘+2) i 


E E |((v'=+'^)dv'f)) 


E , E|((v‘+^)..(^!f’“"))" 


• Wli. 


(40) 

By substituting ([39]) and (HOl) in ([38]) . and equating factors at wu when 
^ = J) we obtain the equality (l37[) . ■ 

By analogy can be proved the following lemma. 

Lemma 2.18 //WA = U = {uij) € with IndU = k and A G M, then 
(AI + U^+2)j, (uf ^) = + ... + 

where = (U^+^)j, (u^-^^) and r‘f^'^ = E ((U^"*"^)! • (u^^^)) « for 

“G/s.nfij 

all s = l,n — 1, i = 1, m, and j = 1 , n. 
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Theorem 2.8 If A G W € with k = ma^{Ind{AW), IndfWA)} 

and rank(AW)^ = r, then the W-weighted Drazin inverse ^ 

(j^mxn respect to W possess the following determinantal representations: 


d,W 
o.o- = 


E , 

lSG.Jr, m 


|((AW)"f 


-fk) 




E 


(AW) 


k+2 d 

h 


(41) 


E |((WA))+2(nf))) 

6 ^r,n{j} _ _ 

E |(WA)^+ 2 g 

<AGlr, n 


(42) 


where is the jth column ofA^ = (AW)^A for all j = 1, ...,m and 
is the ith row of = A(WA)^ for all i = 1,..., n. 


Proof. At first we shall prove (|4T]) . By (IMl) . 


Ad^w = lim ( Aim + (AW) 
A^o ' 


k+2 


-1 


Let 


Aim + (AW) 


k+2 


-1 


1 


(AW)^A. 


/ Lii L 21 

Li2 L 22 


det (Aim + (AW) 


k-\-2 'i 


Lml 

L>m2 


\ I'lm I^2m ■ ■ ■ h^mm J 


where Ljj is a left ij-th. cofactor of a matrix Aim + (AW)^'*'^. Then we have 
(Aim + (AW)^+2)"^ (AW)^A = 


/ 


det(AI^ + (AW)''+2) 


III ( J ■. IIL 

j: E 

S = 1 5 = 1 


4A . -{k) 

E Ls2V, 


(k) 

s2 

(k) 


^ r -ik) \ 
E PslVsn 


5=1 

m 


5=1 


5 I 


V T V T 

/ V ^sm^g\ 2-^ ^sm^g2 


E Ls2V^^ 


S =1 


(fc) 


V s= 


S =1 


E Ls2vfn 

5=1 


AG J -{k) 
E hisruVsn 
s=l 
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By ([M]) . we obtain 

! |(AI™+{AW)'=+2) 

|(AI™+(AW)'=+2)| 

^d,w = lim 

|(AI^+(AW)'=+2)^(v(f)) 
V |(AI^+{AW)'=+2)| 

By Theorem 12.21 we have 
' AI^ + (AW)'=+2 


|(AI^+(AW)'=+2) 

|(AI™+(AW)'=+2)| 


|(AI^+{AW)'=+ 2 )^(vW) 
|(AI™+{AW)'“+2)| / 


(43) 


= A”^ + + dsA™-' + ... + dm, 


where 4 = X) + (AW)^+ 2 ) i 


l3eJs,7 


is a sum of principal minors of 


(AW)^"''^ of order s for all s = l,m — 1 and dm = |(AW)^+^|. 

Since rank(AW)^+^ = rank(AW)^+^ = rank(AW)^ = r, then dm = 
dm-i = ... = dr+i = 0. It follows that det (AI^ + (AW)^+^) = A™ + 
diA™-i + daA™-^ + ... + drX^-^. 

By Lemma [2.171 


Aim + (AW) 


k+2 


_{fc) 


=+ ... + 


for i = l,m and j = l,n, where ^ 


all s = 1 , m — 1 and Cm^ ^ = 


(AW) 


k+2 




((AW) 


k+2 I-(k) 


for 


We shall prove that ^ = 0, when k > r+lfori = l,m and j = 1, n. By 


\k+2 (-(k) 


Tie mm a 12.151 ((AW)*^^^ < r, then the matrix ((AW)^^' 

has no more r linearly independent columns. 

Consider ((AW)*^A2 (vy^^ ^ when /3 G Js,ra{i}. It is a principal sub¬ 
matrix of ((AW)^A2 order s > r -|- 1. Deleting both its i-th row 

and column, we obtain a principal submatrix of order s — 1 of (AW)*^+^. 
We denote it by M. The following cases are possible. 

• Let s = r -|- 1 and det M 7 ^ 0. In this case all columns of M are right- 
linearly independent. The addition of all of them on one coordinate to 
columns of ((AW)^A^ ^ their right-linear independence. 

Hence, they are basis in a matrix ((AW)*^A^ 

column is the right linear combination of its basis columns. From this, 
I ((AW)''+^ ^ ^ ^ Js,n{i} and s = r 1 . 
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If s = r + 1 and detM = 0, than p, [p < r), columns are basis in 


M and 
well. 


in ((AW)^, 


fc+2 1 


■] 


Then 


((AW)"+2 


-Sk) 


= 0 as 


If s > r + 1, then detM = 0 and p, {p < r), columns are ba¬ 


sis in the both matrices M and 


|((AW)'=+2 


-Ak) 


^1=0. 




-Sk) 


Therefore, 


Thus in all cases we have | i'O) ^ ~ when /3 G Js,m{i} and 


r -|- I < s < m. From here if r -|- 1 < s < m, then 




and Cm ^ = det 


Hence, 


((AW)‘+2 (v'f 

(AI + (AW)^+2) . (v 


E |((AW)'-'« 


-Sk) 


= 0 , 


= 0 for z = 1, m and j = 1, n. 


-ik) 


= + ... + for i = 


l,m and j = l,n. By substituting these values in the matrix from 
obtain 


we 


Ad,w = lim 

A—>0 


—1 _j_ —T- 

\”^+di\”^-^+...+dr\"^-^ 

dj- 


w 

df 

(mn) 


41")a"*-i+...+41")a"‘-’' ^ 

\”'+d]_\'^-^ + ...+dr\'^-'^ 

^(mn)yri-lj_ yn-r 

X^+di\”^-^+...+dr\'^-^ 

\ 

/ 




and dr = Y1 (A^+^) ^ 

m 


where ^ 

jd^Jr, m{^} 

Thus, we have obtained the determinantal representation of Ad^w by (l4T]l . 
By analogy can be proved 


3 Cramer’s rules for generalized inverse solutions 
of systems of linear equations 

An obvious consequence of a determinantal representation of the inverse 
matrix by the classical adjoint matrix is the Cramer rule. As we know. 
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Cramer rule gives an explicit expression for the solution of nonsingular linear 
equations. In [35], Robinson gave an elegant proof of Cramer rule which 
aroused great interest in finding determinantal formulas for solutions of some 
restricted linear equations both consistent and nonconsistent. It has been 
widely discussed by Robinson |3S|, Ben-Israel |3S]) Verghese m, Werner 
[38] . Chen [39], Ji [ID] ,Wang [H], Wei [33] , 

In this section we demonstrate that the obtained analogues of the ad¬ 
joint matrix for the generalized inverse matrices enable us to obtain natural 
analogues of Cramer’s rule for generalized inverse solutions of systems of 
linear equations. 


3.1 Cramer’s rule for the least squares solution with the min¬ 
imum norm 


Definition 3.1 Suppose in a complex system of linear equations: 

A • X = y (44) 

the coefficient matrix A G and a column of constants y = (t/i, ..., ym)'^ £ 

C™. The least squares solution with the minimum norm of ^4^ is the vector 
G C"’ satisfying 

llx^ll = min < ||x|| I IIA • x — y|| = min || A • x — y|| >, 

" " xgC" [ xeC" J 

where C” is an n-dimension complex vector space. 

If the equation (I44p has no precision solutions, then x® is its optimal ap¬ 
proximation. 

The following important proposition is well-known. 


Theorem 3.1 ]23f The vector x = A+y is the least squares solution with 

the minimum norm of the system 0- 


Theorem 3.2 The following statements are true for the system of linear 
equations 0- 


i) //rank A = n, then the components of the least squares solution with 
the minimum norm x*^ = (x)*,..., x^ are obtained by the formula 


where f = A*y. 


det(A*A).j (f) 
det A* A 


(V/ = l,n) , 


(45) 
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a) If rank A = r < m < n, then 


Xa = 


E , 

h&Jr,n{j} 


((A*A),,(f)) 


dr (A*A) 


(Vj = l,n) 


(46) 


Proof, i) If rank A = n, then we can represent A"*" by (jlGD . By multi¬ 
plying A+ into y we get (jlSll . 

ii) If rank A = k < m < n, then A"*" can be represented by (dH). By 
multiplying A"*" into y the least squares solution with the minimum norm 
of the linear system (I44p is given by components as in (I46h . ■ 

Using ()13p and (|I7I) . we can obtain another representation of the Cramer 
rule for the least squares solution with the minimum norm of a linear system. 


Theorem 3.3 The following statements are true for a system of linear 
equations written in the form x • A = y. 

i) //rank A = m, then the components of the least squares solution = 
yA+ are obtained by the formula 


0 ^ det(AA*)j. (g) 
* det A A* 


(Vi = 1, m) , 


where g = yA*. 


ii) If rank A = r < n < m, then 


E |((AA*),.(g)): 

dr (AA*) 


(Vi = 1, m) . 


Proof. The proof of this theorem is analogous to that of Theorem 13.21 

Remark 3.1 The obtained formulas of the Cramer rule for the least squares 
solution differ from similar formulas in ]3(k\3^\4l^ . They give a closer ap¬ 
proximation to the Cramer rule for consistent nonsingular systems of linear 
equations. 


3.2 Cramer’s rule for the Drazin inverse solution 

In some situations, however, people pay more attention to the Drazin inverse 
solution of singular linear systems [42H45] . 
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Consider a general system of linear equations (|44l) . where A € and 

X, y are vectors in C”’. i?(A) denotes the range of A and N{A) denotes the 
null space of A. 


The characteristic of the Drazin inverse solution A^y is given in [26] by 
the following theorem. 


Theorem 3.4 Let A € with Ind{A) = k. Then A^y is both the 


unique solution in R{A^) of 


A^+^x = A’^y, 


(47) 


and the unique minimal P-norm least squares solution of 

Remark 3.2 The P-norm is defined as ||x||p = ||P“^x|| for x G C”', where 


P is a nonsingular matrix that transforms A into its Jordan canonical form 


m- 


In other words, the the Drazin inverse solution x = A^y is the unique 
solution of the problem: for a given A and a given vector y G R{A^), find 
a vector x G R{A^) satisfying Ax = y, IndA = k. 


In general, unlike A'^y, the Drazin inverse solution A^y is not a true 


solution of a singular system (I44p . even if the system is consistent. However, 
Theorem 13.41 means that A^y is the unique minimal P-norm least squares 
solution of (1411) . 

The determinantal representation of the P-norm least squares solution 
a system of linear equations ()44p by the determinantal representation of the 
Drazin inverse ()20p is introduced in |46] . 

We obtain Cramer’s rule for the P-norm least squares solution (the 
Drazin inverse solution) of (|44l) in the following theorem. 


Theorem 3.5 Let A G with Lnd{A) = k and ra, uk = rankA^ = 


r. Then the unique minimal P-norm least squares solution x = (xi,..., Xn)"^ 
of the system [44\ ) is given by 



Vi = 1, n, 


(48) 


Xi = 


E (A^+')? 


P^Jr,n 


where f = A^y. 
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Proof. Representing the Drazin inverse by ()26ll and by virtue of Theorem 
we have 


Xl 


/ n \ 

' E disVs ' 


X= I ... 1 =A^y = 


V Xr, 


dr (A*^+1) 


s=l 


E dnsUs 

\s=l / 


Therefore, 


Xi = 


1 ” 

E (A‘«(a<?)) 


4(A^+i)£^ 


/3e Jr.nfj} 


• Vs = 


dr (A^+1) 
1 

dr (A^+1) 


n 

E E (A‘y(a<?)) 


/3GJr,n{i} S = 1 


• Vs = 


n 

E E (Ay‘(a'?.!,.)) 


/ 3 £Jr,n{i} « = 1 

From this (j48p follows immediately. ■ 

If we shall present a system of linear equations as, 

xA = y, 


(49) 


where A G with Ind{A) = k and rankA^"*"^ = rankA^ = r, then by 

using the determinantal representation of the Drazin inverse (|25p we shall 
obtain the following analog of Cramer’s rule for the Drazin inverse solution 
of 


Xi = 


E , (Af+Eg 

OL^Ir, 7i{'i} 

j: |(a‘+i): 

OL^Ir.n 


Vi = 1, n, 


where g = yA^ 


3.3 Cramer’s rule for the W-weighted Drazin inverse solu¬ 
tion 

Consider restricted linear equations 

WAWx = y, (50) 
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where A G W G h = Ind{AW), k 2 = Ind{WA) with 

y G (WA )*^2 and rank(WA )*'2 = rank(AW)*^i = r. 

In [33], Wei has showed that there exists an unique solution A^^wY of 
the linear equations (|50l) and given a Cramer rule for the W-weighted Drazin 
inverse solution of (l5Up by the following theorem. 

Theorem 3.6 Let A, W be the same as in Suppose that U G 

and V* G be matrices whose columns form bases for N{(WA)^‘^) 

and A((AW)^i), respectively. Then the unique W-weighted Drazin inverse 
solution X = (xi, ...jXm of f5dl) satisfies 




WAW 

V 


where i = 1, 2,m. 

Let k = max{A:i, A: 2 }. Denote f = (AW)^A • y. Then by Theorem 12.81 
using the determinantal representation ()4ip of the W-weighted Drazin in¬ 
verse Afi^w^ we evidently obtain the following Cramer’s rule the W-weighted 
Drazin inverse solution of (1501) . 



(51) 


Xi = 


k+2 /3 
h 


m 

where i = 1,2, ..,m. 


E (AW)^+^J 


Remark 3.3 Note that for ^57f) unlike Theorem \3.6[. we do not need auxil¬ 
iary matrices U and V. 


3.4 Examples 

1. Let us consider the system of linear equations. 


2xi — 5x3 + 4x4 = 1) 


7xi — 4x2 — 9x3 + 1-5x4 = 2, 
3xi — 4x2 + 73:3 — 6.5x4 = 3, 


(52) 


Xi — 4x2 + 12X3 — 10.5X4 = 1- 
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/2 0 -5 4 \ 

7 _4 _9 15 

The coefficient matrix of the system is the matrix A = „ , _ ^ ^ 

^ 3-4 7 -6.5 

\1 -4 12 -10.5/ 

We calculate the rank of A which is equal to 3, and we have 


A* = 



7 

3 

1 ^ 



/ 

63 

-44 

-40 

-11.5\ 

0 

-4 

-4 

-4 

A * A 



-44 

48 

-40 

62 

-5 

-9 

7 

12 


L - 


-40 

-40 

299 

-205 

V4 

1.5 - 

-6.5 - 

-10.5yi 



V 

-11.5 

62 

-205 

170.75/ 

we 

obtain entries 

of A+ 

by ([E]): 








63 

-44 

-40 



63 

-44 

-11.5 


d3{A*A) = 

-44 

48 

-40 

+ 


-44 

48 

62 

+ 



-40 

-40 

299 



-11.5 

62 

170.75 




63 

-40 

-11.5 


48 

-40 

62 

+ 

-40 

299 

-205 

+ 

-40 

299 

-205 


-11.5 

-205 

170.75 


62 

-205 

170.75 


= 102060, 


2 

-44 

-40 


2 

-44 

-11.5 


2 

-40 

-11.5 

0 

48 

-40 

+ 

0 

48 

62 

+ 

-5 

299 

-205 

-5 

-40 

299 


4 

62 

170.75 


4 

-205 

170.75 


^11 — 

= 25779, 

and so forth. Continuing in the same way, we get 


A+ = 


1 


102060 


( 25779 -4905 20742 

-3840 -2880 -4800 

28350 -17010 22680 


-5037 \ 
-960 
-5670 

V39558 -18810 26484 -13074/ 


Now we obtain the least squares solution of the system (I52p by the matrix 
method. 


/x?\ 


/ 25779 

-4905 

20742 

-5037 \ 



^0 

^2 

1 

-3840 

-2880 

-4800 

-960 


2 

\xy 

102060 

28350 
\ 39558 

-17010 22680 

-18810 26484 

-5670 

-13074yi 


3 

VJ 


1 

102060 


/ 73158 \ 
-24960 
56700 
\ 68316 / 


/ 12193 \ 
17010 


^71 
9 

5693 

\ 8505 / 


32 






































Next we get the least squares solution with minimum norm of the system 
(I52p by the Cramer rule (j46p . where 


f = 


Thus we have 


tO - 

Xi — 


1 


/2 

0 

-5 

V4 


102060 


7 

-4 


3 

-4 


1 

-4 


/1\ 

2 


/ 26 \ 
-24 


-9 

7 

12 


3 


10 ■ 

1.5 

-6.5 

-10.5/ 

V) 

V- 

-23/ 

26 

-44 

-40 


26 

-44 

-11.5 

-24 

48 

-40 

+ 

-24 

48 

62 

10 

-40 

299 


-23 

62 

170.75 


+ 


- 

^2 — 


X 4 = 


+ 


1 


26 -40 -11.5 

10 299 -205 

23 -205 170.75 


102060 


- 

— 


+ 


1 


102060 


73158 12193 


102060 17010’ 


63 

26 

-40 


63 

26 

-11.5 


-44 

-24 

-40 


-44 

-24 

62 


-40 

10 

299 


-11.5 

-23 

170.75 



-24 -40 62 

10 299 -205 

-23 -205 170.75 


-24960 

102060 


416 

1071’ 


+ 


1 


102060 



48 

-40 

-24 


-40 

299 

10 


62 

-205 

-23 


68316 5693 


102060 8505 


63 

-44 

26 



63 

26 - 

11.5 

-44 

48 

-24 

-h 


-40 

10 - 

-205 

-40 

-40 

10 



-11.5 

-23 170.75 

48 - 

-40 

62 - 

-24 62 

10 -205 

-23 170.75 

N 



56700 

5 





102060 

9’ 


63 

-44 

26 



63 

-40 

26 

-44 

48 

-24 

-h 

-40 

299 

10 

-11.5 

62 

-23 



-11.5 

-205 

-23 


+ 


+ 
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2. Let us consider the following system of linear equations. 


Xi - X2 + X3 + X4 = 1, 

X2-X‘i + X4^ = 2 , 

Xi - X2 + X3 + 2 X 4 = 3 , 

Xi - X2 + X3 + X4 = 1 . 


(53) 


The coefficient matrix of the system is the matrix A 
It is easy to verify the following: 


/I -1 1 1\ 

0 1-11 
1-112 
\1 -1 1 1 / 



-4 

4 

3\ 


/lO 

-14 

14 

10\ 

0 

1 

-1 

0 

, A3 = 

-1 

2 

-2 

-1 

4 

-5 

5 

4 


13 

-18 

18 

13 

V3 

-4 

4 

y 


\io 

-14 

14 

10/ 


and rank A = 3, rank A^ = rank A^ = 2. This implies k = Ind{A) = 2. We 
obtain entries of A^ by ()26l) . 





10 


14 


10 

14 


10 

«2(A j : 


-1 


2 

+ 

13 

18 

+ 

10 


2 


-2 


2 


-1 



18 

13 

+ 

-18 

18 

+ 

-14 

10 

+ 


14 

10 



3 


14 


3 

14 



3 10 

dll — 

0 2 

+ 

4 

18 

+ 


3 10 


10 

10 

= 8 , 


= 4, 


and so forth. 


Continuing in the same way, we get A^ = 


Now we obtain the Drazin inverse solution x 
Cramer rule (148 p . where 


/ 0.5 

0.5 

-0.5 

0.5 \ 

1.75 

2.5 

-2.5 

1.75 

1.25 

1.5 

-1.5 

1.25 

\ 0.5 

0.5 

-0.5 

0.5 / 


of the system (l53l) by the 


g = aV = 

/3 -4 4 3\ 

0 1-10 

4-554 


2 

3 


/10\ 

-1 

13 


CO 

1 

CO 




\io/ 
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Thus we have 



1 / 

10 

-14 



10 

14 


10 

10 

Xl = 

8 ( 

-1 

2 


+ 

13 

18 

+ 

10 

10 


1 / 

10 

10 




-1 

-2 


-1 

-1 

X2 = 


-1 

-1 

+ 


13 

18 

+ 

10 

10 


1 / 

' 10 

10 




2 

-1 


13 

13 

X3 = 

8 ( 

13 

13 

+ 


-18 

13 

+ 

10 

10 



10 

10 




2 

-1 


18 

13 

X4 = 

8 ( 

10 

10 

+ 


-14 

10 

+ 

14 

10 


_ 1 

“ 2 ’ 

= 1 , 

= 1 , 

_ 1 
“ 2 ' 


4 Cramer’s rule of the generalized inverse solu¬ 
tions of some matrix equations 

Matrix equation is one of the important study fields of linear algebra. Linear 
matrix equations, such as 

AX = C, (54) 

XB = D, (55) 

and 

AXB = D, (56) 

play an important role in linear system theory therefore a large number of 
papers have presented several methods for solving these matrix equations 
[IMI]. In [52], Khatri and Mitra studied the Hermitian solutions to the 
matrix equations (1571) and (|56p over the complex field and the system of the 
equations (|57|) and (l5^ . Wang, in [53l[5l|, and Li and Wu, in [55| studied 
the bisymmetric, symmetric and skew-antisymmetric least squares solution 
to this system over the quaternion skew field. Extreme ranks of real matrices 
in least squares solution of the equation (f56ll was investigated in [56| over 
the complex held and in m over the quaternion skew held. 

As we know, the Cramer rule gives an explicit expression for the solution 
of nonsingular linear equations. The Robinson’s result ( [35|) aroused great 
interest in hnding determinantal representations of least squares solution as 
analogs of the Cramer rule for the matrix equations (for example, |58H60] i. 
The Cramer rule for solutions of the restricted matrix equations ()57[) . (I55I) 
and (f56]l was established in [6TH63] . 

In this section, we obtain analogs of the Cramer rule for generalized 
inverse solutions of the matrix equations (f57|) . (f55]l and (f63l) without any 
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restriction, namely for the minimum norm least squares solution and Drazin 
inverse solution. 

We shall show numerical examples to illustrate the main results as well. 

4.1 Cramer’s rule for the minimum norm least squares solu¬ 
tion of some matrix equations 

Definition 4.1 Consider a matrix equation 


(57) 


AX = B, 


where A G G given, X G is unknown. Suppose 

Si = {X|X G ||AX - B|| = min}. 


Then matrices X G such that X G S'! are called least squares solutions 

of the matrix equation ( f57l ). IfKis = ™iiixeSi ll^ll? then ^ls is called the 
minimum norm least squares solution of ( |57| ). 

If the equation ([57)1 has no precision solutions, then X.ls is its optimal 
approximation. 

The following important proposition is well-known. 

Lemma 4.1 ( The least squares solutions of (5^ are 


X = A+B + (I„ - A+A)C, 


where A G C™'’^”',B G are given, and C G is an arbitrary 

matrix. The least squares minimum norm solution is 'X.^s = A’^^B. 

We denote A*B =: B = (kj) G C’"’'". 

Theorem 4.1 (i) If rank A = r < m < n, then for the minimum norm 

least squares solution X^s = (xij) G for all i = l,n, j = 1, s we 
have 



(58) 


/S^Jr, n 

(ii) If rank A = n, then for all i = l,n, j = 1, s we have 


_E (A-A)g 



(59) 


det(A*A) 


where b j is the jth column of B for all j = 1, s. 
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Proof. i) If rank A = r < m < n, then by Theorem 12.31 we can represent 
the matrix A”*" by (|lip . Therefore, we obtain for all i = l,n, j = 1, s 


„ E ((A*A).,(ay)? 

E l n{*} 

k=i E (A*A) ^ 


iSeJr. ' 


■ ^kj — 


E Er=i ((A*A).aa*,)) 

l3eJr,n{i} 


■b 


kj 


E (A*A) ^ 

/S^Jr, n 


Since E ^*k^kj = 


( E ^Ik^kj \ 
E 

k 


\ 


= bo', then it follows 


E ^nk^kj 

k 

(ii) The proof of this case is similarly to that of (i) by using Corollary 12.11 


/ 


Definition 4.2 Consider a matrix equation 


XA = B, (60) 

where A G G are given, X G is unknown. Suppose 

S 2 = {X| X G ||XA - B|| = min}. 

Then matrices X G such that X G S '2 are called least squares solutions 

of the matrix equation [WD . IflK^s = ™iiixeS 2 II^IL then is called the 
minimum norm least squares solution of (6^. 

The following lemma can be obtained by analogy to Lemma l4.ll 

Lemma 4.2 The least squares solutions of IgQj) are 

X = BA+ + C(I^-AA+), 

where A G G are given, and C G is an arbitrary 

matrix. The minimum norm least squares solution is'X.^s = BA'*'. 

We denote BA* =: B = (kj) G 
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Theorem 4.2 (i) If rank A = r < n < m, then for the minimum norm 

least squares solution 'Kls = {xij) € for all i = l,s, j = l,m 

we have 


Xij — 


E 

a€lr,m{j} 


AA*) . (b,. 


E 


(AA*)g| 


(61) 


OtGir, m 


(a) If rank A = m, then for all i = l,s, j = l,m we have 

det(AA*)j. (bi.) 
det(AA*) 


(62) 


where bj. is the ith row of B for all i = 1, s. 


Proof, (i) If rank A = r < n < m, then by Theorem 12.31 we can represent 
the matrix A"*" by (I12p . Therefore, for all i = 1, s, j = 1, m we obtain 


n n 

Xij = bik • 

k=l k=l 


E . 


((AA*)^. (a^.)) 


a 

a 


E l(AA*)g 

m 


ELi 


'ik 


E . 


AA* 



a 

a 


E l(AA*)g 

m 


Since for alH = 1, s 


Y^b,,al = (Eb^kaU Ebikal, ••• E = b,„ 

k 

then it follows m- 

(ii) The proof of this case is similarly to that of (i) by using Corollary 12.11 


Definition 4.3 Consider a matrix equation 

AXB = D, (63) 

where A G G (f^mxq given, X G is unknown. 

Suppose 

Ss = {X| X G II AXB - D|| = min}. 
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Then matrices X € such that X G S's are called least squares solutions 

of the matrix equation ^6^. IfKis = minxeSsll^ll? then is called the 
minimum norm least squares solution of i63\) . 


The following important proposition is well-known. 

Lemma 4.3 ( fMjl) The least squares solutions of [6^) are 

X = A+DB+ -h (I„ - A+A)V W(Ip - BB+), 

where A G G G are given, and {V, W} C are 

arbitrary quaternion matrices. The minimum norm least squares solution is 

Xi5 = A+DB+. 

We denote D = A*DB*. 


Theorem 4.3 (i) If rank A = ri < n and rankB = r 2 < p, then for the 

minimum norm least squares solution 'Kls = {xij) G of l[6M) we 

have 


j h 


E (A-A),(dT* 

E (A*A)^ E |(BB*) 


/ 3 GJr 




(64) 


or 


E (BB*)^.,(df) “ 
0‘&Ir2,p{j} 

E (A*A)^ E l(BB* 


(65) 


/SeJi-j,' 


(XGir 


where 


d" = 


d^ = 


E 

.aG/r2,p{i} 


:bb*), (di. 


(BB*)^, (d„. 

«e/r2,p{i} 


H T 


( 66 ) 


E (A*A),,(d,i) ^ (A*A)., (d.,^ ^ 

7i\i^ 


(67) 

are the column-vector and the row-vector, respectively, dj. is the ith 
row of Y) for all i = l,n, anddj is the jth column of T) for all j = l,p. 
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(a) //rankA = n and rankB = p, then for the least squares solution 
'^LS = {xij) G of we have for all i = l,n, j = l,p, 


Xi 4 — 


det((A*A).i (dB 


or 


where 


d? - 


df := 


det(A*A) • det(BB*) ’ 

det((BB*),. (d^)) 
det(A* A) • det(BB*) ’ 

det((BB*),-. (di.)) ,...,det((BB*),-. (d,. 
det ((A*A).i fd.i)) ,... ,det ((A*A).j (d. 


1 T 


( 68 ) 

(69) 

(70) 

(71) 


are respectively the column-vector and the row-vector. 


(Hi) /frank A = n and rankB = r 2 < p, then for the least squares solution 
^LS = (xij) £ of we have 


(72) 


or 


det ( (A*A). j ( d^ 


det(A*A) ^ |(BB*); 

(y.^Ir‘2 ,p 

E (BB*)^. (d.A) 
06 /^ 2 ,pit} 


det(A*A) E l(BB^ 


a€l, 


\a I ? 

)a I 


(73) 


r2,p 


where d® is [UBi) and df- is ([7l| ). 


(iiii) /frank A = ri < m and rankB = p, then for the least squares solution 
^LS = (xij) G of (E3) we have 

det((BB*),. (d^)) 


Xij — 


or 


Xij — 


E (A*A)P •det(BB*; 

,n 

E (A-A),.(d = )? 

/3G Jri,n{*} 


(74) 


E (A*A)^ det(BB*) 

,n 


(75) 


where d® is fTOl ) and df- is ([67l ). 
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Proof. (i) If A € B € and n < n, r 2 < p, then by Theorem 12.31 

the Moore-Penrose inverses A"'' = (afj\ € and B+ = (ti) € 


possess the following determinantal representations respectively, 


E (A*A) ^ 


d y /3 


“«= 


E (A-A) 

, n 




E |(BB*),.(bE; 

aelr2,p{j} 

E l(BB*)g| 

Q:G/ro.T) 


(76) 


Since by Theorem 14.31 Xrg = A+DB^, then an entry of ^ls = ( 


Xij) IS 


Xij — 


i2if2^tkdks]btj. 


(77) 


S = 1 \fc = l 


Denote by d.^ the sth column of A*D =: D = (dij) € for all s = l,q. 

It follows from E ^*kdks = d.^ that 


E (A*A) , (a* ^ ^ 

m m ^ ' z .1 \ . 


^tk^ks— 

k=l 


, 71 Id 


KkJ 13 


k=l 


E (A*A) ^ 


E E (A*A),aa 

n{i} ^ — 1 


A y^.k) p 


/3eJri , n 

dks 


dks — 


E , (A*A)^d 


E (A*A) 

/3eJri , n 


E (A*A) 

/3eJri , n 


(78) 


Suppose Os. and e. g are respectively the unit row-vector and the unit column- 
vector whose components are 0, except the sth components, which are 1. 
Substituting (ITHI) and (|76|) in ([771) . we obtain 


Xij — 


E 

S = 1 


E (A*A),Jd, 

/3GJri.n{d 


E |(BB*),,(b,7): 

"GTj.pDI 


E (A*A) ^ 

jB^Jr-^ , n 


E l(BB* 

CH^lr 


r2^P 
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Since 


ds = '^e,idis, h* ='^b*tet., '^disb*^ = du, 


(79) 


i=i 


t=i 


S = 1 


then we have 


Q P n 

EEE E (A*A).aeO^« E |(BB*),-.(eO; 

S = 1 t=l 1 = 1 l3^Jr^,n{i} aG/r2,p{j} 


0 GJri , n 


E (A*A)^ E l(BB*) 


/3 


aG/i 


T2>P 


p n 

E E 

t=l«=l/3eJri,n{*} 


EE E (A*A).ae,,)? rftt E l(BB-),.{e,): 




E (A*A) 2 E KBB^* 


Denote by 




dt ■= 


(80) 


(A-A).,(d,^^ 


E E (A-A),ae,)? 

/ = 1 /3GJri,n{j} 


dit 


the tth component of a row-vector = (d^,..., d^) for all t = l,p. Sub 
stituting it in (ISOll . we have 


Xjd — 


E4 E l(BB*),,(eO 

t=l a^Ir2,p{j} 


l3GJri , 7 


E (A*A)g E l(BB 


O'G.Ir 


Since E d^^t. = d^, then it follows 
t=i 

If we denote by 


dfj--=Y.dit E l(BB*),,(eO“|= E | (BB*),-. (d,) (81) 

t=l a£lr 2 ,p{j} ae^r 2 ,p{i} 
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the /th component of a column-vector d^- = (d®,for all I = l,n 
and substitute it in ()80p . we obtain 


Xij — 


E E 

I 1 ^SJri, n'iO 


(A*A).ae,)^ 


d 


B 

b 


E 

, n 


{A*A) ^ 


E l(BB*)« 

Q^G/r2 ,p 


Since ^ = d^-, then it follows ([Mil . 

1=1 

(ii) If rank A = n and rankB = p, then by Corollarv l2.1l A+ = (A*A)~^ A* 
and B+ = B* (BB*)~^. Therefore, we obtain 




^LS 

= (A*A)~iA*DB 

* 

(BB* 

)■' = 



Xi2 ... 

Xlp ^ 




fLti 

■ 



2^21 

X22 ... 

X2p 


1 


Li2 

Lt2 • 

• ^n2 






det(A* A) 





V 

^nl 

Xn2 • • • 

Xnp ) 





Lt • 

rA , 

• ^nn/ 


/ dll 

di2 ■ ■ ■ 

dim \ 


/R^i 

RB .. 



d2i 

d22 ■ ■ ■ 

d2m 


1 


Ru 

pB 

rX22 . . 

• R% 






det(BB*) 






\dnl 

dn2 . . 

dnm 

/ 


\R% 

R\ •• 

dB , 

. ^pp/ 


where dij is ijth entry of the matrix D, is the ijth cofactor of (A* A) for 
all i,j = l,n and is the ijth cofactor of (BB*) for all i,j = l,p. This 
implies 

E ^ki ( E dksR^s ) 

X - ■ — _E ( 82 ) 

det(A*A) • det(BB*) ’ 

for all i = 1, n, j = l,p. We obtain the sum in parentheses and denote it as 
follows 

p 

Y^dksRfs = det(BB*),-. (d^.) := df^, 

S = 1 

where d^, is the A:th row-vector of D for all k = l,n. Suppose d^- := 
... ,d^j^ is the column-vector for all j = l,p. Reducing the sum 

n _ _ 

we obtain an analog of Cramer’s rule for (I63p by (I68p . 

k=l 
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Interchanging the order of summation in ([8^ . we have 

p / n ^ \ 

E [ELtidkARfs 

_ S=1 \k=l _ / 

det(A*A) • det(BB*)' 

We obtain the sum in parentheses and denote it as follows 

n 

Y^Ltdks = det(A*A),, (d.,) =: 4, 

k=l 


where d.^ is the sth column-vector of D for all s = l,p. Suppose d^ := 


4 5 • • •) df- ) is the row-vector for alH = 1, n. Reducing the sum ^ 


A 
i . 

A dB 


s=l 


we obtain another analog of Cramer’s rule for the least squares solutions of 
(pjl by (IMl)- 

(iii) If A € B € and ri = n, r 2 < p, then by Remark 


IQ and Theorem 12.81 the Moore-Penrose inverses 

B+ = 


A+ = 




€ C’’ 


and 


respectively, 


€ C‘^^P possess the following determinantal representations 

, det(A*A).da*' 


% = 


det (A*A) 




E |(BB*),.(bE 

aG/r2,p{j} 


E 


BB* 


(83) 


a€l, 




Since by Theorem 14.31 Xrg = A+DB^", then an entry of = (x^) is 

(177)1 . Denote by d.^ the sth column of A*D =: D = (dij) G for all 

s = l,g. It follows from E^*kdks = d.^ that 

k 


^tkdks — ^ 


det (A*A) , (a*,) ^ det (A*A) , (^d.„^ 


k=l 

Substituting 


k=l 

and 


det (A*A) det (A*A) 

in (1771) . and using (1791) we have 


S = 1 


det (A*A) 


E l(BB*) 

(y^Ir2,p 
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Q P n 

E E Edet(A*A).Je.;H.6:t E l(BB*),.(eO 

S It 1/ — 1 ,p 0} 

det(A*A) E l(BB*)g| 


ael. 


^2’P 


p n 


E Edet(A*A) ,(e,0 du E |(BB*),.(eO 

t=ll = l 0‘£lr2,p{j} 

det(A*A) E l(BB*)g| 

Oi^Ir^ ,p 

If we substitute (j81|) in (185 1) . then we get 


(85) 


Xij — 


E det(A*A) . (ej) df- 
1=1 _ _ 

det(A*A) E l(BB*)“ 

OL^Iro , 7 ) 


Since again E ~ d^-, then it follows (f72]l . where d^- 


IS 


1=1 

If we denote by 


dA 


det (A* A). i (d. t j =Y det (A* A) Jej) du 
1=1 1=1 

the tth component of a row-vector = {df[, ...,d^) for all t = l,p and 
substitute it in (f85|) . we obtain 


E4 E |(BB*),.(e02l 

t=i ae/rj.plil 


det (A* A) E l(BB 


OL \ 


Oi^lr 


Since again E d^^t. = d^, then it follows (f73l) . where d^ is (f7T]l . 

t=i 

(iiii) The proof is similar to the proof of (iii). ■ 


4.2 Cramer’s rule of the Drazin inverse solutions of some 
matrix equations 

Consider a matrix equation 

AX = B, (86) 

where A € with IndA = A:, B G ([^nxm given and X € -g 

unknown. 
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Theorem 4.4 ( Theorem 1) If the range space R(B) C R{A^), then 
the matrix equation with constrain RfX.) C R(A^) has a unique solution 

X = A^B. 


We denote =: B = {kj) G 

Theorem 4.5 If rankA^+^ = rank A^ = r < n for A G then for 

Drazin inverse solution X = A^B = (xjj) G of ^8^) we have for all 

i = '^,n, j = l,m, 


Xij — 


E 

/ 3 eJr.n{*} 



E 

0&Jt, n 


(A‘+i) ^ 


/3 

h 


(87) 


Proof. By Theorem 12.71 we can represent the matrix A^ by (1261) . There¬ 
fore, we obtain for all i = l,n, j = 1, m, 


S=1 S=1 


E , 

n{0 


fc + 1 (^{^)\\ ^ 


E 

n 


■ bsj — 


E . EEi 

j3SJr, n{T 


fc+1 I 

h 


JSJ 


^^Jr, n 


Since ^ sS^^bgj = 


e4‘>.j 


E O-niKj 


= h j, then it follows ([87)1 . 


Consider a matrix equation 


XA = B, (88) 

where A G with IndA = /c, B G are given and X G is 

unknown. 
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Theorem 4.6 ( f6^ , Theorem 2) If the null spaee A^(B) D N{A^), then the 
matrix equation [8^) with constrain NfK) D N{A^) has a unique solution 

X = BA^. 


We denote BA^ =: B = (kj) G 

Theorem 4.7 //rankA^+^ = rankA^ = r < m for A G then for 

Drazin inverse solution X = BA^ = (xij) G of I88\) . we have for all 

i = j = l,m, 


E 

aelr,m{j} 


(b- 


E l(A'=+') 

m 


(89) 


Proof. By Theorem 12.71 we can represent the matrix A^ by (1251) . Therefore, 
for all i = l,n, j = 1, m, we obtain 


Y1 = Y1 


^ ^is^sj — 
s=l s=l 


E 

Ol&Ir,m{j} 


, k+1 ( (k) 


J ■ 


a 


E l(A‘+‘) 

OL^Ir, m 


Er=i hk E 

06 T.m{i} 


k+1 f^{k) 
j ■ 


E l(A‘+‘) 

Oi^Ir, m 


Since for alH = 1, n 


hisaif^ = ( E bisa^si E hsa 


Ak) 

^s2 


E bisdirri^ = bj , 


then it follows 

Consider a matrix equation 


AXB = D, (90) 

where A G with IndA = fci, B G with IndB = k 2 and 

D G are given, and X G is unknown. 

Theorem 4.8 ( Theorem 3) //i?(D) C R{A^^) and N(D) D 
k = max{ki,k 2 }, then the matrix equation with constrain ii(X) C 

R{A^) and A(X) D A(B*=) has a unique solution 

X = A^DB^. 
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We denote =: D = (d^) 

Theorem 4.9 //rankA^i+^ = rankA^i = n < n for A G and 

rankB*^2+i = rankB^^ = r 2 < m for B G i^m^m^ then for the Drazin 
inverse solution X = A^DB^ =: {xij) G of ^ we have 


Xjd — 


E 


d-?- ? 


deJri ,n 


E (a'»+‘)2 E |(b‘ 


fc2 + l'l 


OL^lr 


(91) 


or 


where 


Xij — 


E 

(y.^Ir2 ,m {j} 


b; 


fc2 + l 


(di^) 


deJi-i. 


E {A‘>+‘)^ 


E |(B‘=+') 


OtGir 


(92) 


= 


E 

OL^Ir2 ,m {j} 


Bfc2+i (d- ' " 


E 

Q^Ir2 ,m {j} 


B 


fc2 + l 


H T 


d 1 “ 

On. I a 


(93) 


= 


E iAE‘(d.)^|,.... E 


fci+i 


d 


_/3eJri.n{*} OG^ri.n{*} 

are the column-vector and the row-vector, dj, and d ^- are respectively the 
ith row and the jth column o/D for all i = l,n, j = l,m. 


Proof. By 


B^ = 


and ()25p the Drazin inverses A^ = 


.D 


G C" 


and 


G possess the following determinantal representations, 


respectively, 


E 


fci+1 /3 


AE." 


/3 


^ij 


D _ {j} 

hij — 


E (A‘>+1)^ 

/3^Jr^ , n 


E l(B^^+')gl 


(94) 


OtSlr 
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Then an entry of the Drazin inverse solution X = A^DB^ =: {xij) G 
is 


m / n 


Xij — 


s=l \t=l / 

, fcjv fS — 


(95) 


Denote by d.^ the sth column of A^D =: D = (dij) € for all s = l,m. 

It follows from ^ Si^dts = d. s that 


t=i t=i 


E 

Ti{i} 


A':**' ^ 


E (A‘>+‘) ^ 

jS^Jr-^ , n 


df.ei - 


E E 




• dts E 

pG.Jri, n{'i} 


(d 


E (A^1+^) B 


13 & Jr I , n 


E (A^i+I) ^ 


/3^Jr-^ 

1 


Substituting (IM|) and (fM)l in ([95]), we obtain 


(96) 


Xij — 


E 

S=1 


E 




Q;G/r2 ,m 0)' 


E (A'-''+‘) ? 


0GJr-i , 


E l(B"^+') 


CX-^Iv 


Suppose e^, and e. ^ are respectively the unit row-vector and the unit column- 
vector whose components are 0, except the sth components, which are 1. 
Since 


d.6 = '^e^idis, J^disbft^'’ = du, 




1=1 


t=l 


S=1 


then we have 


Xij — 


m m n 


E EE E 

S = li=l; = l/3GJri.n{i} 


fel+1 


e.z 


disOgt 


OS/rj {j} 


E 


E (A"i+')^ E l(B"^+') 

/3G(/r]^,n OL^Ii 


r2,Tn 
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(97) 


m n 

EE E 

A‘y‘(e,)^ 

dit E 




OL^Ir2 







E |(B* 2 + 1 ) 


Denote by 


di--= 


E . 

, n{^} 


A-'i+l f3,) ? 


E E !A‘r‘(e,)? 

1=1 0^Jr^,7i{i} 


dit 


the tth component of a row-vector = (d^,..., for all t 

Substituting it in (IWl) . we obtain 


= 1, m. 


Xnj - 


E4 E , 
^ 1 




E E l(B"=+') 


ael, 


r2,m. 


Since ^ = df", then it follows ([92]). 

t=i 

If we denote by 


t—l aG/rj.mfi} 


E |B-=+'(dz.): 

fiG/7‘2 {j} 


the lih component of a column-vector d® = (d®,...,d^)^ for all I 
and substitute it in (fW]) . we obtain 


= 1, n 


Xin - 


E E 

I 1 /3G<7r]^, n{0 


, ^1 + 1 


(e.O 


% 


l3eJ, 


E E l(B^=+') 




aei, 


r2,m 


n 

Since Yl ^-idy = d^-, then it follows ([OT]1 . ■ 


4.3 Examples 

In this subsection, we give an example to illustrate results obtained in the 
section. 
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1. Let us consider the matrix equation 

AXB = D 


( 98 ) 


where 


A = 


/I i i\ 


/I z 1\ 

i -1 -1 

, B = ( ‘ ' , D = 

i 0 1 

0 10 
V-1 0 -i) 

V-1 * 1 / 

1 i 0 

\0 1 ij 


Since rank A = 2 and rankB = 1, then we have the case (ii) of Theorem 
Ol We shall find the least squares solution of ([Ml) by ([M|) . Then we have 


A*A = I -2i 3 2^ , BB* = , D = A*DB* = -1 


and |(BB*) “I = 3 + 3 = 6, 
2 


E (A* A)? 

0 &J 2 , 3 


= ^"*'-2z 3j+^"M2 3)+^"M-3 z ? 


= 10 . 


By ([70]), we can get 

/ 1 2z 3z\ 

Since (A*A) ^ (d®) = j —i 3 2 1, then finally we obtain 

\-i 2 3/ 


xii = 


E (A-A)i(dB)^ 


/3G J2,3 

Similarly, 


E (A-A)g E l(BB> 


drti_i, 3)+drt(j;, ^ 

60 60' 


«e/i,2 


, , —i 2i\ . , (—i 3i 

detl j 3j+det(_3 3 

60 




Xu = 
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det 


X21 = 


—2i —i 


+ det 


-i 2 
-i 3 


60 


det 


X22 = 


3 -i 
-2i -1 


+ det 


-1 2 
-1 3 


60 


det 


X31 = 


—3i —i 


+ det 


3 -i 
2 -i 


60 


det 


X32 = 


3 -i 
-3i -1 


+ det 


3 -1 
2 -1 


60 


2i 

2 

i 

1 


2. Let us consider the matrix equation (I98p . where 

, B = 

^1 i Oy 

We shall find the Drazin inverse solution of (jOSp by (I64p . We obtain 
A2 = 




4 

0 

0\ 

/ 8 

0 

2 - 2z 

0 

0 

, A^ = 4 - 4z 

0 

-2 - 2i 

0 

0/ 

\-4 - 4i 

0 

1 


—z 

i 3 — i\ 


B^ = 


1 

-1 l + 3z . 



^—3 + i 3 — i 3 i / 

Since rank A = 2 and rankA^ = rankA^ = 1, then ki = IndA = 
ri = 1. Since rankB = rankB^ = 2, then k 2 = IndB = 1 and r 2 = 2. 
we have 


and X; « 

/3eJi,3 


-4 

D = A^DB = ( -2 + 2i 
2 + 2i 

= 8 + 0 + 0 = 8 , 


4 

2 - 2i 


4- 4i 


-2-2f -4-4zy 


E |(B^)S| = 

, 3 


det 


1 -1 


+ det 


3 — z 3 + 1J 


—i 3 — i 
—3 + i 3 + z 


0 + (-9 - 9z) + (9 - 9z) = -18z. 


2 and 
Then 
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By ([66]) . we can get 


'12 - 12i^ 
= ( -12i 

-12 




^-12 + 12i\ / 8 

12i , = -12 - 12i I . 

12 / \-12 + 12i. 


'12 - 12i 0 

Since (d^) = j —12i 0 0 | , then finally we obtain 

-12 0 0 > 


Xu = 


E |a’. (d=) I 

/36Ji,3{1} 


12-12f 1 + i 


E (A?t E l(B=):i 8^ (-IS) 12 


/3gJi- 


CtGl2,3 


Similarly, 


Xl2 = 


-12 + 12i -1-i 


8 ■ (-18i) 


12 ’ 


3)13 = 


8-(-18z) 18 


-12f 1 


12f 


8-(-18z) 12 


1 -12-12f 1-i 

, 3)23 = 


12 


3)31 = 


Then 


8 • (-18i) 


i -12 

, 3)32 = 


8 • (-18z) 12 ’ 

-12 + 12i -1-i 


12 


X = 

is the Drazin inverse solution of 


8.(- 

■18f) 

12- - 

l-l-i 

— 1—i 


12 

12 

18 \ 

1 

1 

1-i 

12 

12 

12 . 

i 

i 

-1-z / 

12 

12 

12 / 


8 • (-18i) 


12 


5 An application of the determinantal representa¬ 
tions of the Drazin inverse to some differential 
matrix equations 

In this section we demonstrate an application of the determinantal repre¬ 
sentations of the Drazin inverse to following differential matrix equations, 
X' -|- AX = B and X' -I- XA = B, where the matrix A is singular. 
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Consider the matrix differential equation 

X' + AX = B 


( 99 ) 


where A G B € i^rixu given, X G i^rixu unknown. It’s well- 

known that the general solution of (|99p is found to be 


X(t) = exp 


-At 


exp^‘ dt ) B 


If A is invertible, then 


exp^* dt = A ^ exp^* -|-G, 


where G is an arbitrary n x n matrix. If A is singular, then the following 
theorem gives an answer. 


Theorem 5.1 ( Theorem 1) If A has index k, then 
I exp^* dt = A^ exp^* +(I-AA^)t 


A A^ 9 A*’ ^ i, 

^+y‘+iir‘ +■■■ + —* 


+G. 


Using Theorem 15.11 and the power series expansion of exp we get an 
explicit form for a general solution of (|99l) 

X(t) = 

{a^ + (I - AA^)t (l - ft + + g} B. 

If we put G = 0, then we obtain the following partial solution of (1991) . 


X(t) = A^B + (B - A^AB)t - i(AB - A^A2B)t2 + ... 


(- 1 )'' 


(A^-^B - A"A^B)t' 


k\ 


D \ki 


( 100 ) 


Denote A^B =: B^') = (bf^) G for all I = 1, 2k. 

Theorem 5.2 The partial solution I1100\} . X(t) = (xij), possess the follow- 
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ing determinantal representation, 


Xoj - 


E 

n{0 




i fc+i 




+ I bij — 


E 

0€.Jr, n{*} 




0GJr, n ' 

, E 

1 I “^(l) 0^Jr, n{i} 

“2 1 % ~ 


E (A'=+1)^ 


..+i(g(.+2))) 


3eJr, 

(-l)fe I PeJr,n{i} 


E (A'“+1)^ 


t^ + 


Afe + l/'g(2fe)'l'| /3| 


’)) 


fc! 




E (A'=+1) 

0^Jr, n 


( 101 ) 


for all i,j = 1, n. 


Proof. Using the determinantal representation of the identity A^A (I32p . we 
obtain the following determinantal representation of the matrix A^A”^B := 


yij = = Y1 

S=1 t=l 0&Jr,n{i} 


E, 

/3Gt/r,n{^} 


E 

S = 1 


Ay (a.(‘+‘))) 


E“ 

t=l 


St ^tj 


E (A-'+i)? 

P^Jr,n 


n 

E 

t=i 


): 

• btj E 

(Ay (b‘^”>)) 1 

E 

(A^+i)' 

E 

n 

(A*+i) 



for all i^j = 1, n and m = 1, /c. From this and the determinantal representa- 
tion of the Drazin inverse solution (|87p and the identity (I32p it follows (llOip . 


Corollary 5.1 If IndA = 1, then the partial solution of [9^) . 

X(t) = (xij) = A^B + (B - A9AB)t, 
possess the following determinantal representation 


Xij — 


E 

0^Jr, nfd 




/ 


I^SJr, n 

for all i,j = 1, n. 


E (a'^)2 


+ 


bij 


E . 

f^GJr, n{^} 


A2, b 


:(2) 


E (A2)^ 

0^Jr, • 


( 102 ) 
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Consider the matrix differential equation 

X' + XA = B 


( 103 ) 


where A G B G are given, X G is unknown. The general 

solution of (llOSp is found to be 


X(t) 



If A is singular, then an explicit form for a general solution of (11031) is 


X(t) = 

B IA^ + (I - AA^)t (l - 41 + + g} • 

If we put G = 0, then we obtain the following partial solution of (11031) . 


X(t) = BA^ + (B - BAA^)t - i(BA - BA‘^A^)t‘^ + ... 
(BA^-^ - BA'=A^)t^. 


(104) 


Denote BA^ =: B^^^ = (bfj) G for all I = 1,2A:. Using the determi- 

nantal representation of the Drazin inverse solution (I89p . the group inverse 
(I30p and the identity (|31l) we evidently obtain the following theorem. 

Theorem 5.3 The partial solution 1^104^ , X(t) = (xtj), possess the follow¬ 
ing determinantal representation, 


E 

{j} 




E |(A'=+I)g| 


+ 1 bij — 


E 


k+1 f 


)): 


E |(A'=+I)g| 


-- I - 


E 

Oi^Ir,n{j} 


k + 1 ( •L(fc+2)') 


E |(A'=+I)g 

i k + 1 


e + ... 


. / ^ 

( — 1)^ I 0‘^lr,n{j} 

~^r~ 1 




E |(A'“+1); 


for all i,j = 1, n. 


Corollary 5.2 If IndA = 1, then the partial solution of I1103\} . 

X(t) = {xij) = BA^ + (B - BAA^)t, 
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possess the following determinantal representation 


Xij — 


E 


A] [h) 


( 1 ) 


/ 


E l(A2) 

(AGlr,n 


+ 


bij 


E 

a€lr,n{j} 


A2 (b} 


( 2 ) 


E l(A2) 

OtGir.n 


for all i,j = 1, n. 


5.1 Example 

1. Let us consider the differential matrix equation 

X' + AX = B, 

where 


(105) 


1 -1 


'1 i V 


A = 


—i i , B = i 0 1 


.-1 1 


.1 i 0 , 


Since rank A = rankA^ = 2, then k = IndA = 1 and r = 2. The matrix 
A is the group inverse. We shall find the partial solution of (jin5l] by (I102p . 
We have 


A2 = 


I 

i 

3-i\ 

2-i 

2i 

0\ 

1 

-1 

l + 3i 

, b(^) = AB = 1 + 2i 

-2 0 

V— 3 + i 3 — i 

3 + i/ 

V 1 + i 

i 

0/ 


^2 - 2z 2 + 3f 

B(2) = A^B = ( 2 + 2i -3 + 2f 0 | . 

T+5i -2 Oy 


and 


E (A^) 

3 


2^ h 
P 


det 


1 -1 


+ det 


-1 1 + 3A 

3 — i 3 + w 


+ det 


—i 3 — i 
—3 + i 3 + i 


0 + (-9 - 9f) + (9 - 9f) = -18z. 

/ 2 — i i 3 — i\ 

Since (A^) ^ = j 1 + —1 1 + 3i j and 

\l + i 3-i 3 + i / 

/2 — 2i i 3 — i\ 

(A2) ^ = 2 + 2f -1 l + 3z 

\l + 5f 3-i 3 + i/ 


57 
















then finally we obtain 


E |(A^(bW))^| / 

E 1 (a^(£(=))) ^1 


/5e32,3{l}'^ ^ 


E (A2); 

^e-^2,3 




E (A2); 

P£J2,3 


t = 


3-3t 

-18i 


+ ( 1 -^]t = 


-18i 


6 • 


Xii = 

Similarly, 

—3 + 3f / 9 + 9A —1 —i 1 + i . , 

“^■2 = ^18-+ •-^1*=-^ + —*. x.3 = 0 + (l-0)f = t, 


3 + 3f / -18 \ -1 + i 

-18i V ) 6 ’ 

—3 — 3z / —9 + 9A 1 — i 1 + i 

X22 = -TTT^ + ( 0-7^ ) t = — -^- —t, X23 = 0 + (1 - 0) t = t, 


-18f 


-18f 


6 


-12i ( -m\ 2 

3^31 = —777 +1-rwr ] t = -, 


9 + 3i 


-18i 

-18 

-18i 


t = 


-18z 

-1 + 3i 
6 


3^33 = 0 + (0 - 0) t = 0. 


Then 


I 1 + i —1 —i + (3 + 3f)t t 
X = -( —1 + i 1 — z + (3 + 3z)t t 
4 -l + 3z 0 


6 


is the partial solution of (11051) . 


References 

[1] E. H. Moore, On the reciprocal of the general algebraic matrix. Bulletin 
of the American Mathematical Society 26(9) (1920) 394”395. 

[2] A. Bjerhammar, Application of calculus of matrices to method of least 
squares; with special references to geodetic calculations, Trans. Roy. Inst. 
Tech. Stockholm 49. (1951). 

[3] R. Penrose, A generalized inverse for matrices. Proc. Camb. Philos. Soc. 
51 (1955) 406-413. 


58 
























[4] M. P.Drazin, Pseudo-inverses in associative rings and semigroups, The 
American Mathematical Monthly 65(7) (1958) 506-514. 

[5] K.M. Prasad, R.B. Bapat, A note of the Khatri inverse, Sankhya: Indian 
J. Stat. 54 (1992) 291-295. 

[6] I. Kyrchei, Analogs of the adjoint matrix for generalized inverses and 
corresponding Cramer rules, Lin. Multilin. Alg. 56(4) (2008) 453-469. 

[7] I. Kyrchei, Analogs of Cramer’s rule for the minimum norm least squares 
solutions of some matrix equations, Appl. Math. Comput. 218 (2012) 
6375-6384. 

[8] I. Kyrchei, Explicit formulas for determinantal representations of the 
Drazin inverse solutions of some matrix and differential matrix equations, 
Appl. Math. Comput. 219 (2013) 1576-1589. 

[9] I. Kyrchei, The theory of the column and row determinants in a quater¬ 
nion linear algebra. Advances in Mathematics Research 15, pp. 301-359. 
Nova Sci. Pubh, New York, 2012. 

[10] I. Kyrchei, Cramer’s rule for quaternion systems of linear equations, J. 
Math. Sci. 155 (6) (2008) 839-858. 

[11] I. Kyrchei, Cramer’s rule for some quaternion matrix equations, Appl. 
Math. Comput. 217(5) (2010) 2024-2030. 

[12] I. Kyrchei, Determinantal representation of the Moore-Penrose inverse 
matrix over the quaternion skew held, J. Math. Sci. 180(1) (2012) 23-33. 

[13] I. Kyrchei, Determinantal representations of the Moore-Penrose inverse 
over the quaternion skew held and corresponding Cramer’s rules, Lin. 
Multilin. Alg. 59(4) (2011) 413-431. 

[14] 1. Kyrchei, Explicit representation formulas for the minimum norm least 
squares solutions of some quaternion matrix equations, Linear Algebra 
Appl. 438(1) (2013) 136-152. 

[15] 1. Kyrchei, Determinantal representations of the Drazin inverse over the 
quaternion skew held with applications to some matrix equations, Appl. 
Math. Comput. 238 (2014) 193-207. 

[16] 1. Kyrchei, Determinantal representations of the W-weighted Drazin 
inverse over the quaternion skew held, Appl. Math. Comput. 264 (2015) 
453-465. 


59 



[17] X. Liu, G. Zhu, G. Zhou, Y. Yu, An analog of the adjugate matrix for 

( 2 ) 

the outer inverse Mathematical Problems in Engineering, Volume 

2012, Article ID 591256, 14 pages, doi:10.1155/2012/591256. 

[18] R. B. Bapat, K. P. S. Bhaskara, K. M. Prasad, Generalized inverses 
over integral domains. Linear Algebra Appl. 140 (1990) 181-196. 

[19] A. Ben-Israel, Generalized inverses of marices; a perspective of the work 
of Penrose, Math. Proc. Gamb. Phil. Soc. 100 (1986) 401-425. 

[20] R. Gabriel, Das verallgemeinerte inverse eineer matrix, deren elemente 
einem beliebigen Korper angehoren, J.Reine angew math. 234 (1967) 
107-122. 

[21] P. S. Stanimirovic’, General determinantal representation of pseudoin¬ 
verses of matrices. Mat. Vesnik 48 (1996) 1-9. 

[22] A. Ben-Israel, On matrices of index zero or one, SIAM J. Appl. Math. 
17 (1969) 1118-1121. 

[23] R. A. Horn, C. R. Johnson, Matrix analysis. Gambridge etc., Cambridge 
University Press, 1985. 

[24] P. Lancaster, M. Tismenitsky, Theory of matrices, Acad. Press., New 
York 1969. 

[25] C.F. Van Loan, Generalizing the singular value decomposition, SIAM 
J. Numer. Anal. 13 (1976) 76-83. 

[26] Y. Wei, H. Wu, The representation and approximation for the weighted 
MoorePenrose inverse, Appl. Math. Comput. 121 (2001) 17-28. 

[27] Stephen L. Campbell, Carl D. Meyer Jr., Weak Drazin inverses, Linear 
Algebra Appl. 20(2) (1978) 167-178. 

[28] P.S. Stanimirovic’, D.S. Djordjevic’, Full-rank and determinantal repre¬ 
sentation of the Drazin inverse, Linear Algebra Appl. 311 (2000) 131-151. 

[29] Carl D. Meyer Jr. , Limits and the index of a square matrix, SIAM J. 
Appl. Math. 26(3) (1974) 506-515. 

[30] R. E. Cline, T. N. E. Greville, A Drazin inverse for rectangular matrices, 
Linear Algebra Appl. 29 (1980) 53-62. 


60 



[31] Y. Wei, Integral representation of the W-weighted Drazin inverse, Appl. 
Math. Comput. 144 (2003) 3-10. 

[32] Y. Wei, C.-W. Woo, T. Lei, A note on the perturbation of the W- 
weighted Drazin inverse, Appl. Math. Comput.149 (2004) 423-430. 

[33] Y. Wei, A Characterization for the W -weighted Drazin inverse and a 
Cramer rule for the W - weighted Drazin inverse solution, Appl. Math. 
Comput. 125 (2002) 303-310. 

[34] Z. Al-Zhour, A. Kiligman, M. H. Abu Hassa, New representations for 
weighted Drazin inverse of matrices, Int. Journal of Math. Analysis 1(15) 
(2007) 697-708. 

[35] S.M. Robinson, A short proof of Cramer’s rule. Math. Mag. 43 (1970) 
94-95. 

[36] A. Ben-Israel, A Cramer rule for least-norm solutions of consistent lin¬ 
ear equations. Linear Algebra Appl. 43 (1982) 223-226. 

[37] G.C. Verghese, A Cramer rule for least-norm least-square-error solution 
of inconsistent linear equations. Linear Algebra Appl. 48 (1982) 315-316. 

[38] H. J. Werner, On extension of Cramer’s rule for solutions of restricted 
linear systems, Lin. Multilin. Alg. 15 (1984) 319-330. 

[39] Y. Chen, A Cramer rule for solution of the general restricted linear 
equation, Lin. Multilin. Alg. 34 (1993) 177-186 

[40] J. Ji, Explicit expressions of the generalized inverses and condensed 
Cramer rules. Linear Algebra Appl. 404 (2005) 183-192. 

[41] G. Wang, A Gramer rule for minimum-norm (T) least-squares (S) so¬ 
lution of inconsistent linear equations. Linear Algebra Appl. 74 (1986) 
213-218. 

[42] G. Wang, A Gramer rule for finding the solution of a class of singular 
equations. Linear Algebra Appl. 116 (1989) 27-34. 

[43] G. Chen, X. Chen, A new splitting for singular linear system and Drazin 
inverse. J. East China Norm. Univ. Natur. Sci. Ed., 3 (1996) 12-18. 

[44] A. Sidi, A unified approach to Krylov subspace methods for the Drazin- 
inverse solution of singular nonsymmetric linear systems. Linear Algebra 
Appl. 298 (1999) 99-113. 


61 



[45] Y. M. Wei, H. B. Wu, Additional results on index splittings for Drazin 
inverse solutions of singular linear systems, The Electronic Journal of 
Linear Algebra 8 (2001) 83-93. 

[46] P.S. Stanimirovic. A representation of the minimal P-norm solution. 
Novi Sad J. Math. 3(1) (2000) 177-183. 

[47] H. Dai, On the symmetric solution of linear matrix equation, Linear 
Algebra Appl. 131 (1990) 1-7. 

[48] M. Dehghan, M. Hajarian, The reflexive and anti-reflexive solutions of a 
linear matrix equation and systems of matrix equations, Rocky Mountain 
J. Math. 40 (2010) 825-848. 

[49] F.J. Henk Don, On the symmetric solutions of a linear matrix equation. 
Linear Algebra Appl. 93 (1987) 1-7. 

[50] Z.Y. Peng and X.Y. Hu, The generalized reflexive solutions of the ma¬ 
trix equations AX = D and AXB = D, Numer. Math. 5 (2003) 94-98. 

[51] W.J. Vetter, Vector structures and solutions of linear matrix equations. 
Linear Algebra Appl. 9 (1975) 181-188. 

[52] C.G. Khatri, S.K. Mitra, Hermitian and nonnegative definite solutions 
of linear matrix equations, SIAM J. Appl. Math. 31 (1976) 578-585. 

[53] Q.W. Wang, Bisymmetric and centrosymmetric solutions to systems 
of real quaternion matrix equations, Comput. Math. Appl. 49 (2005) 
641-650. 

[54] Q.W. Wang, The general solution to a system of real quaternion matrix 
equations, Comput. Math. Appl. 49 (56) (2005) 665-675. 

[55] Y. Li, W/ Wu, Symmetric and skew-antisymmetric solutions to systems 
of real quaternion matrix equations, Comput. Math. Appl. 55 (2008) 
11421147. 

[56] Y. H. Liu, Ranks of least squares solutions of the matrix equation 
AXB=C, Comput. Math. Appl. 55 (2008) 1270-1278. 

[57] Q.W. Wang, S. W. Yu, Extreme ranks of real matrices in solution 
of the quaternion matrix equation AXB=C with applications, Algebra 
Colloquium, 17 (2) (2010) 345-360. 


62 



[58] G. Wang, S. Qiao, Solving constrained matrix equations and Cramer 
rule, Appl. Math. Comput. 159 (2004) 333-340. 

[59] Y. Qiu, A. Wang, Least-squares solutions to the equations AX = B, XC 
= D with some constraints, Appl. Math. Comput. 204 (2008) 872-880. 

[60] Y. Yuan, Least-squares solutions to the matrix equations AX = B and 
XC = D, Appl. Math. Comput. 216 (2010) 3120-3125. 

[61] G. Wang, Z. Xu, Solving a kind of restricted matrix equations and 
Cramer rule Appl. Math. Comput. 162 (2005) 329-338. 

[62] C. Gu, G. Wang, Condensed Cramer rule for solving restricted matrix 
equations, Appl. Math. Comput. 183 (2006) 301-306 

[63] G.J. Song, Q.W. Wang, H.X. Chang, Cramer rule for the unique solu¬ 
tion of restricted matrix equations over the quaternion skew field, Com¬ 
put. Math. Appl. 61 (2011) 1576-1589. 

[64] Z.Xu, G. Wang, On extensions of Cramer’s rule for solutions of re¬ 
stricted matrix equtions. Journal of Lanzhou University 42 (3) (2006) 
96-100. 

[65] S. L. Campbell, C. D. Meyer, JR. and N. J. Rose, Applications of the 
Drazin inverse to linear systems of differential equations with singular 
constant coefficients, SIAM J. Appl. Math. 31 (1976) 411-425. 


63 



